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ABSTRACT 


Thermophotovoltaic cells are a good candidate for use in high efficiency radioiso- 
tope thermoelectric generator (RTG) power devices for deep space missions. This thesis 
examines the use of Silvaco Virtual Wafer Fabrication Software as a tool for designing 
and optimizing TPV cells for different possible spectra. It gives results for GaSb and In- 
GaAs cells optimized to the AMO spectrum which closely match published data as well 
as hypothetical cells optimized to the spectrum of a 1300K blackbody. 


THIS PAGE INTENTIONALLY LEFT BLANK 


vi 


= 


Il. 


TABLE OF CONTENTS 


INTRODUCTION wisivessscscscedesecunssvadivisdesacecsssessdadstedsnesuatestsdescbeenseussbeudsdaddeaveitnsadsceseves 1 
SOLAR CEES isc sicaisise <ics wean’ sueahvanisvsbussuedsssanvadssveussaasssenscanssovassonssudssvgavscdacbonseounvenss ey 3 
A. FATS POR Y spssiscsessiss ceasnsies sastsaevcscvssssed sone tvecsnyned dusssauveconvveved soneddessnened casstunvsseusvensien 3 
B. THE PEN JUNCTION 6 sscsiinsessscasdsniscssscesssousssdasevssssacsssesiestevexssteasssdsssuseessoonsaniene 4 
Cc PHOTOVOLTAIC CEL DU vissscessssicossnctsacensessbisssinsssassevasooussnacesdeseboopousssbunsexaers 8 
D. ENCIDEN PE SPHC TRUM viscvssicstedcovetcevedsasdsvsnidcdevtedestersosvtsovsccsdeccsssistcsoedsteateck 8 
E. POWER OU TP Pisses sssssnssvecesussunssopececsssguvssassoursecstvenassesedvecsdescessdsvensdesesesss 11 
F, EFFICIENCY IMPROVEMENTS .........ccsccssssssssscssscssssescseccsscssscsessoessesssess 13 
1. ROPE CDG Hh osaccancecdicesucsenicacicicaWineasbcahaeonscskas Seacaoeesdaanisasedaanesceeunsenes 13 

2. Low Bmeray PROMS ss ccisssicesescisateatestssecvsiasseisevasatvacusecincesadisvessactaisivrs 14 

3. High Knerey- PHO tO oes és icecesiecesvssnesccqusstvsens scostedesiedcovedsseceessseecevesedeors 14 

4. COIL ROSISEANCE sevcspenceseedecscuicssaskacsdastedssiiustovscdsucceuttvckaveddstedssesatinrinassoees 14 

5. Material Derects sii iitesiscssacccivitersdccsdesssioutevicssescdessveseyscvesdevcobontersasesads 14 

6. SSID AATTIG sec dceess escadarcatagsacacdensesexeaea tes desneas as teleueaseemaaeanieeten eae race 15 
THERMOPHOTOVOLTAIC CELLS .........ccsccssssssssscssscesssoscsessssssesssossesssscsesseesses 17 
A. PETS POV iivsscscvontcvis doa Gacapsasubscdecuducovsnd dviadstonsosssosdsvedestienapsibeviadodetetessescatedsvades 17 
B. INCIDENT SPECTRUM 'ssssscessssvassonssvecnsenes sessssucsdontvenas sna svavadepevcnssionsseusdocsss 18 
1. Blackbod y¥: Radiata Qi ccssenssassseavsescdenves soseuvsssiccsion dssssedvdesscecnacdstedvdaesssees 18 

Zs Radioisotope Thermoelectric Generator Spectrum...............sceeee 21 

Cc, THERMOPHOTOVOLTAIC SYSTEM. ..........scssscssssssssccssscsssssssscssescsesessnes 22 
1. PRCA SOUP sssseticisyudssadscadsecvnituvsndvtelevcionsdeidisescaceseusadiedebedeveseiccstadeisiae 23 

2; DEE CH VE HMMMELED ss, sudesscascscaesssseusadcesdacsosasteesdesstedessesesvuacssdaccaissturisesseeee 23 

3. BUG sevsivsssevshconciscatsebs eudsoadscesviduvs ddstosevbvontsvacdsescassswebavaddeteducssentusudivadds 23 

4. Photo- Electric Converter .............scccsccscsssssccscscssccsssscsscssscscsssssscseseeees 25 
RADIOISIOTOPE THERMOELECTRIC GENERATOR ...........cccsscsssssssssseseeees 29 
A. PETS POR gs cssecsioctsiscsnscades cesabsedeseesdovetcavsddauddadssscecvideesceddusutbovscdesddceasdchcasecevsevs 29 
B. BUED SUPPLY. sssvcsssssssssissacusccsesoesssonscucssvenstcasaseucslccuavenascensdvesssssstedsesedisescesense 32 
C. SABE DY . ccvssevas sana cuicenscavdseassauvssonsvenes sonedsavasensdcasshecdsdesdnensteassdvsongestecbsssvsceusssssss 34 
1. ACCIGONE: HISKOPY secscssainesadsareuaisctiacsescas cctadantasdaasaarssdscbaisesbedeascaceahacsaveas 34 

2 Safety: MO ASU OS. 5s ssiiccasiiesiisscuessadaicavesbevavcassercexsisaaiousdeinenaveoesinanetcence 35 

3. OAL OEY SONTAG 5c saedaseacecsoacesvecevesanansdendsxsceesdachedesvatverscsecdcecnneezaicsoszets 36 

a. PTGS, pcan cud ciales csliavisstusd sands asisvcad cudecsas tuxadaseviduahie condvidawedeidss cates 36 

b. BI GST cscsctcssansicesasusstensakels asad caitisea hans voeven ves esas sone Beene 36 

Cc. FOCI 3s cisesa shat sess bcoz nceasd cas anestealasptaav aeundceast sia 36 

d. POT TD GCL a i isgnciiceatinia autem 37 

e. DIMM CTSIONAN. WUC i. desins deceit seein eek courses ehaateed 37 

f STEED csuscsscestnaseessingnnsldesksy sean sivas sskay os dsyanas'vueshon wae puncnsdeanon ooetes 37 

g. TAP SO FLA RTIIONIS vvcnciscisssasicuitiivatedcaas salencadviavased cotlesxcctaetacscoides 37 

D. FUTURE [IMPROVEMENTS scsicsasccssesasvcssscestscnssscssscnsostssbousscessvensecossuntsocnse 38 
1. Dynamic Isotope Power Systems ...........ssccsssccssssccssssccssssccssssscssssscees 38 


vil 


2; Alkaline Metal Thermal-to-Electric Conversion ...........ccccccccssssssees 40 


3. Thermophotovoltaic Cells ..............ccscsscccssssssccssccsscccssssescccesssssscssoeees 41 

V. SIMULA TION SOP TP WARE: yccccisivasscsssncasccabstoessciesecessevcsseaasngeitccbeeueststssvoncebeanvbassou 43 
A. SEIS V ACO si seacas bscbvicdeisvesisscecusssnisatieassa eusedapecvsiosabeockscswseuansusGoncguiversiereenbscbnie 43 

B. DEEEENIN GA CRU sisssioccctsacesssveceseanstscshensscvscstuotsenscescobsaasosesnesttcassckonsbubeseoes 45 

1. Defining Constants iiicssccsscenstescacssscnscenscesnvesscnssevenacsvensnecesceavenonxsvncsonse 46 

2. Defining Mesh Bound aries..............csccccsssccssssccssssccsssscesssscsssscsssescees 46 

3. Building the X- and Y- MeShS ..............cssssssccssssssccsssssccssssssesscesees 47 

4. Defining ROS OMS sccsiscvevscsnctevscveccevacseseoastvbcenvens nevuncenyeuetevsnanotecpeeevecewee 48 

5. PRACing PICCIrOdeS eyes scisecidissantstaiessciidececicaacdiedzucdentdanadenelsiacsvecstecanees 50 

6. Defining Doping Levells............sccccsssscssssssssssscssssssesssssescscssscsssosssssees 50 

a Defining Material Properties ...............ccssccsssscsssssccsssscesssscessescesseccess 50 

8. FrCIGENE S PO CEP UI 5 sakscecdicccsiechecoasuscseqsuvereaesoudacsestuogceousesuccesneeseeveroeneve 51 

9, SOM VLIND scaceu ches denscvesecavavoasetesdestevnssiurunssereecesadwevelaccessuaesevssenanesasererounetevens 52 

10 SOUDCE COGE siseiccinesiedateciatabiseasecsiestecebivcacsentacocetatosescdyecossutinansvessosnesene 53 

c. IEA TTAB scssctistccssiecpocisiesacsdsstabosacsstonctanessessbansocsdualeboceneneniceonesSasidivenstassiawneurs 53 

1 NUGLAROMS saci lesocstiicsciacaccecktsssseacescunecacsecscsussoesusgtsnedssescesvosunsnensdenseavecece 53 

2. PRUE D siss scdcassisdaasivas scceivis Sdasesuvsdaccisuad ssssucdsdeesesddsiisatiuncticlerssbeavectsetseds 53 

VE... (MODEES. AND RESUILCUES wisiesavecatvivessscusesesspsssvedssoussccspetaes cushsesksebetvicnsstavvestboanstenonss 55 
A. MODEL CONETRMA TION sisisccsssisccosssenssosacoossaadentenaibaassdesscoceccsooedsonesasenss 55 

if Callitiny A FSOTIE: caccissocivcicesedesesiceescagcncodeuaveetneteuxbccesedvsnencinesceieuvessvaee 55 

2; Gallium Antimonide ..............scccssscssssscssssscsssssssssssssesscssesscsssssssssessssees 58 

3. Indium Gallium Arsenide .............cccssccssssscssssscssssccsssscsssssssssssscssssseees 60 

B. THERMOPHOTOVOLTAIC CELLS .............sccssssccssssccssssccssscsssssssssssssees 62 

1. Gallic A MGMONIGG  sescsisccisiseciaiccoeticcaciescdesenseseibosesencesenansceceesooaenivees 62 

Zs Indium Gallium Arsenide .............ccccsccssssscssssscssssccssssscssssccssssssscsscees 64 

VII. CONCLUSIONS AND RECOMMENDATIONS. ..........ccccsscssssssssessscsscscssssssssescseee 69 
APPENDEX A? ATLAS CODE scssciciscesissnscveessnvisstensshcavtauteuvecssnecvekcesesnueassouchcaniwasteneevebsciest 71 
A. GALLIUM ARSENEDE: scssscasedscessienossessbsecteancescadaeacssucacessinanssnasbudsasoucdebipecacss 71 

iF DV CAIDVE: sacs soci sstsassoueibenupcussscsoasubonns sotsabteasseaasacannacsoacsnssnpostsatbasibaasnes 71 

2. FreGuencCy RES POnSe vissdiss censsainicesscansescdisdissecsocasiaasabsdacentenssdeteiehionnes 74 

B. GALLIUM ANTIMONIDE............cccccssccsscssssssoscsnsscessccessssescsnecsesscoessscosssesces 80 

1, TeV GQUDVG:sicviciyewvessscocescuosueswchtsueksesyanckaspacssesciusbbecsbhivesseiuae ccapbeneiassons 80 

2. Frequency Response sssiisisSvccstitivcecinseiccusetsscdstaceiiscneceis scindatscenetashininss 84 

C. INDIUM GALLIUM ARSENIDE 10.0.0... ccccsssccccccsccccecscscccsccsescccseesssesensees 87 

1, PAV, GAL VO seatibivvsteincacsbucsestsssuvssuinseuhcasesWesas teas seiesecahvasoxeaboanbsdessbeaustasoous 87 

2. Frequency: RESPODSE sisccisiscsnsciasticacsacetectindssacunscasunedaatesscebsenssennssbbasnace 92 

APPENDIX Bi MATLAB CODE csscisscccsssssvessssnsoncscdocvsesssvensassnvevedseondovtssbovvsoessvensectbsdunbsoonss 97 
A. DEL ASA RUIN GM a, sscdscacaneascccehicieutascabancbingceskantabinesababeccanentuacsbnssoncaneavicwiaaiaaess 97 

B. BUDE WV sIVAs,  cctvca vist hteas so esetodabehcovecusvansuctsncunpsuecssubesshscussotsaasoucsncappcnavatessbabasvess 99 

C. IBV INIAAXCP, Wil, scctvcstactestestsotucasacetasuncesganssucissasnesseausass sacessovsasboutsncevves teal iononbacsbse 99 

D. BEA CK BODY. Mi issiscsiccscessitevsideceevnssstiass cdssscasessssavsedestadessenceosddsessdesisbaviedoses 100 

E. PRE ORESP IME scasisives sctwsssacasveunsasuiuececepveuies snip seu cceucwasdsbegesecan yank unipwehcesbwouebanny 102 


1x 


THIS PAGE INTENTIONALLY LEFT BLANK 


Figure 1. 
Figure 2. 
Figure 3. 
Figure 4. 
Figure 5. 
Figure 6. 
Figure 7. 
Figure 8. 
Figure 9. 


Figure 10. 
Figure 11. 
Figure 12. 
Figure 13. 
Figure 14. 
Figure 15. 
Figure 16. 
Figure 17. 
Figure 18. 
Figure 19. 
Figure 20. 
Figure 21. 
Figure 22. 
Figure 23. 
Figure 24. 
Figure 25. 
Figure 26. 
Figure 27. 
Figure 28. 
Figure 29. 
Figure 30. 
Figure 31. 
Figure 32. 
Figure 33. 
Figure 34. 
Figure 35. 
Figure 36. 
Figure 37. 
Figure 38. 
Figure 39. 
Figure 40. 
Figure 41. 


LIST OF FIGURES 


Poe Tie UTICA ON 5 ah ce ud eas cloe cae cao a ney Seances ase eo a eo ee es 4 
Energy Band Model of a p-n Junction [From Ref. 2.] oo... ee eee eeseeeeeeeeeeeeees 5 
Built-In Potential [After Ref. 2.] .......cccc cc ccccsssssccececececeesessseeececeesesessnseaeeeeees 5 
p-n Junction Under Bias [After Ref. 2.].... ee eee eeeeeeseecseceseeeeeeeeneecnaeeneeneees 7 
Ssinplé SolarCell [Atter Ret: 341 cnnincieive coud tained ected 8 
AMO Spectrum vs. Wavelength [Data for Plot is From Ref. 4.] .... ee 9 
AMO Spectrum vs. Photon Energy [Data for Plot is From Ref. 4.].......0..0.. 10 
Current in a Solar Cell [After Ref. 5]... cccccccesscccecececsesssseceeeeeeeeeesenssaees 11 
POV WIV ies cate estde aa cata eas ess areb eases ap eee aad antes ean nse 12 
Anti Reflective Coating and Textured Cell Surface [From Ref. 6.] ................ 13 
AMO vs. 5800K Blackbody Spectrum ...........eccceeeseeseceecesececeeeeeesteeeeseerenaeeees 20 
L300K: Black Dod y Spur wii sock joes seis ceduczssGhcvniqesigdtsdedeces tee Toceanssged avtuaeieeseedbans 22 
TEP NV SS YSCEM xe ce scarce ates stiaycgenegedeeaseaealoanste sued ivi dgebsaetavan take kbos eeeaiedeharocstedeastans 22 
Resonant Filter [After Ref. 11.) oo. c cc cccccececececececesecesesesesesesueueueeeess 24 
Resonant Filter Frequency Response [From Ref. 11.].... cece eeeeeeeneeeeeeeees 24 
RTG Module [From Ref. 16.] ........cccccccccccccccccccccececececececececececeseseseseseceseeeeeseeess 32 
Pu-235.Power Output [Prom Ref. 17.) sszssecscesenceeaicd es dea deceaeeegiay naetecteter aes ans 33 
“President Shows Atom Generator” [From Ref. 17.] .......ceccccccceeeeeeeeeeees 34 
RIG Safety Measures. [From Ref. 19.) 5.212. de costes iscsconioqustiecceeaieeanivess 35 
Simplified Stirling Engine [From Ref. 20.] ...... cece eeseeeseceseceseeeeeecnaeenseeeeees 38 
110 Watt Electric Stirling Generator [From Ref. 21.] ....... eee eeecceeeseeeeeteeeeeee 40 
AMTEC Schematic [From Ref. 22.] ........cccccccccccccccccccecececececececesesececeseeueueeeeeess 40 
Simulator Block Diagram [From Ref. 23.] ......ceccccesscecescceeeeeeeeeneeeeeeeeeeneeeenes 45 
CP o1 CAL IMIG SIN esses Po excices Se acod cod vracts rs wasee sounatey ba yaacane oy saatucaseccanacaveeanesdencosaparedeteesas 48 
MSH O IONS ticctattenrsiins iaease lead cecsatancsascaseteaeat sancshcotiaecycaunccaseseucie ance as teuee axe 49 
‘GaSb.opt’ [Data for Plot is From Ref. 21]... eee eeececeeceeesseeeeeneeeenneeeenes 51 
GraNS: SOlAT CGM ace Gus sce de uaa tntey vcuie dorccdecaanamadeseusedvetseceayennvensuaeencee seat 56 
RFA eH MRSS LG ics oatss let scab tcslccy cua ceva dee sae searsei esos eu sata moons baad eee nt aecas tae 56 
Spectral Response Of GaAs vecsisescsisvessicccsseasancadeseseessusscdensegsactesuthaeteleodeaeastoedane 57 
RAS DOLD CGE eogretit ck acatdscomcaasvsasndeascamumeseatua scans ssuaecy teacca Mapueieuatincsameneeeatee 58 
GaSb Solar Cell-V C baracteristies jo.5i..0iapiaendaescesidevaveasbasaanelphivedsdeacaheeseeneess 59 
Spectral Response: of Gasp ieiicssesiascissavecaisaseeesasvannadasaceeeeassavaatesdeenecedeaceane eavedde 60 
InGaAs Solar Cell iia asecnsyaceiiasetiaaseeataapseae a aocied OaeiheasdaiiebeaMavaneladsateuens 61 
InGaAs Solar Cell J-V Characteristics .........ceccccesccecseeeceseeecsseceesseeeeseeeenaeeees 61 
Spectral Response Of InGaAs a. .3:1ccicsjcceds sageaasees sdanessussanencs doaatacuaeacsvansacavsseseanee 62 
ASD DP Y Cen cn taacccalsccanha verti tingod aigeh lain sas acganstaes sd tae laninalanednusenieees 63 
GaSb TPV Cell (EV Characteristics i... csisacceidsvsrsadesaccetaesavaatasensaleaspaccavedeaeeade 63 
GaSb Photo@eneration Rate .i.s..scsisassscedaseseaetiaveratedsocsaceanncoaate sobsncdessaveesardenens 64 
Mas TA eI arse ude ieee pecns dag na ease vga nc ous a aoe os 65 
InGaAs TPV Cell J-V Characteristics ..........cecccecssececsscecesececeseceesteeeeseeeenaeeees 65 
InGaAs Photoseneralion Rate ascii. aay loreacsaaaesia ne eisdgchletieaiee aaa 66 


xi 


THIS PAGE INTENTIONALLY LEFT BLANK 


Xil 


Table 1. 
Table 2. 
Table 3. 


LIST OF TABLES 


Température. vs: Peak Wavelengths :2x.nec0 ied ce oe eae et 21 
Maximum Frequency Response of Various Matertals.............ceeeeeseeeseeeeeeees 26 
RTG Missions [After Ref. 15.) .....ccccccccccccccccccecsesesssecececeeecsessnseceeeeeesesesensaaees 31 


xiii 


THIS PAGE INTENTIONALLY LEFT BLANK 


XIV 


ACKNOWLEDGMENTS 


I would like to thank Professor Sherif Michael Professor Karunasiri for their guid- 
ance on this thesis. 


XV 


THIS PAGE INTENTIONALLY LEFT BLANK 


XVI 


EXECUTIVE SUMMARY 


A major limitation in space exploration is the availability of power. While satel- 
lites stationed around Earth can produce an adequate supply of power using solar arrays, 
satellites farther from the Sun require unreasonably massive arrays to produce very little 
power. In order to accommodate the need for power on deep space missions, it is neces- 
sary for satellites to carry a source capable of providing adequate power for the life of the 
mission. This is currently done using radioisotope thermoelectric generators (RTGs). 
Current RTGs consist of a heat source, Pu-238, and a thermocouple that inefficiently 
converts heat into electricity. 

In order to improve the overall efficiency of RTGs, a better thermoelectric con- 
verter is needed to replace the thermocouple. This thesis examines the possible use of 
thermophotovoltaic (TPV) cells in RTGs. Two cells were developed and optimized for 
the spectrum from a 1300K blackbody which simulates the spectrum of an RTG heat 
source. 

Current TPV cells are optimized primarily us ing a triakand-error method. Cells 
are built and tested with different thicknesses and doping levels in order to the most effi- 
cient design. This method is very expensive and time consuming. This thesis presents a 
model that can predict the output of a cell under various spectra. The model can easily be 
changed to simulate different cell thicknesses and doping levels. Additionally, with a 
comprehensive set of material parameters, the model can be changed to simulate cells of 
different materials. 

The models presented in this thesis were built using the Silvaco Virtual Wafer 
Fabrication software package. This software package has been shown by other research- 
ers as a capable tool for modeling advanced solar cell designs, but there is no record of 
any TPV models ever being designed. 

In order to prove the capabilities of this software package, an initial cell was de- 
signed using a very well documented material, gallium arsenide (GaAs). Voltage-current 
characteristics and frequency response data were recorded from this model and compared 
with experimental data from a similar cell. Once the model was verified, more exotic 
materials could be examined. Gallium antimonide (GaSb) and indium gallium arsenide 


XVii 


(InGaAs) were examined for TPV cells because these are two of the most researched ma- 
terials and, as such, there is a significant amount of experimental data to compare with. 

As with the GaAs cell model, the GaSb and InGaAs models were first optimized 
for the AMO spectrum and compared to experimental results. Again, the voltage-current 
relationship and frequency response of the cells were measured and found to be very 
close to available experimental data. Once these models were determined to be accurate, 
they could be optimized for a different spectrum. 

Programs were written using Matlab to run iterations of Silvaco model changing 
layer thicknesses and doping levels in the cells. Other programs gathered data and calcu- 
lated and recorded open-circuit voltage, short circuit current, maximum power, fill factor, 
and efficiency of each cell in order to determine the most efficient design. These results 


were also compared to experimental data on similar cells and found to be very close. 
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I. INTRODUCTION 


As space-based- missions take both manned and unmanned vehicles farther from 
Earth, the need for a better source of power becomes more important. The intensity of 
light from the sun decreases with the square of the distance from it making solar cells an 
inadequate source of power on distant missions because the area of cells required would 
be tremendous. As an example, a 1 m’ solar panel at earth would have to be 25 m? at 
Jupiter and almost 2000 m’ at Pluto to yield the same output power. Batteries, another 
common source of power on satellites are used primarily only as storage devices because 
they are only capable of providing a few hundred watt-hr/kg. This means that the 
amount of batteries required to provide primary power for a long mission would be too 
massive to even launch into space. Current deep space satellites use an alternative source 
of power known as a radioisotope thermoelectric generator or (RTG). 

An RTG consists of a heat source and a means to convert that heat into electricity. 
Existing RTGs accomplish this using thermocouples. Although these are very simple and 
reliable thermoelectric converters, they operate at only a few percent efficiency. In order 
to improve overall efficiency of an RTG, the thermocouple will have to be replaced by a 
more capable thermoelectric converter. 

The most attractive replacement for the thermocouple is currently the thermo- 
photovoltaic (TPV) cell. Research has been ongoing in this field since the 1950’s, but the 
exotic materials necessary for high efficiency cells has only been recently available. 
Here, several models will be developed and optimized for both the sun’s spectrum, AMO, 
and the spectrum from a 1300K blackbody which represents well the spectrum from an 
RTG heat source. 

Due to the fact that TPV cells work on the same principles as solar cells, these are 
discussed in Chapter II. TPV cells are discussed in Chapter III. In Chapter IV, the entire 
RTG module is discussed including the safety precautions necessary when using a radio- 
active heat source. Chapter V discusses the software packages used in these models and 
Chapter VI gives the optimized cell characteristics along with results, and Chapter VII 
contains conclusions and recommendations for future work. ATLAS code used in this 


thesis is in Appendix A and Matlab code is in Appendix B. 
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Il. SOLAR CELLS 


This chapter contains the basic principles of solar cells. This information is im- 
portant to this thesis because thermophotovoltaic cells work on the same principles as 


solar cells. 


A. HISTORY 


The history of the solar cell, also known as a photovoltaic cell, can be traced all 
the way back to 1839 when the then nineteen-year-old French physicist Antoine-César 
Becquerel first discovered the photovoltaic effect [Ref. 1]. While experimenting with 
electrolytic cells containing two different metal electrodes in a solution, he found that 
some combinations of metal and solution would produce a small current when exposed to 


light. 


The first actual solar cell was built in 1877 by Charles Fritts [Ref. 1]. He applied 
a very thin gold film to a piece of selenium. Although this cell operated with less than 
one percent efficiency, it was able to produce a measurable amount of power [Ref. 1]. A 
similar cell was constructed in 1927 consisting of a copper film on piece of copper oxide 


[Ref. 1]. This cell also operated with less than 1% efficiency [Ref. 1]. 


The next breakthrough came in 1941 when Russell Ohl built the first silicon solar 
cell [Ref. 4]. In 1954, Pearson, Chapin, and Fuller expanded on this design and were able 
to build a silicon cell with an efficiency of 6% [Ref. 1]. Since that time, there has been a 
great deal of progress in solar cell design. Gallium arsenide cell designs have been able 
to achieve efficiencies of over 20% and multijunction cells have been built with up to 


37% efficiency [Ref. 1]. 


B. THE P-N JUNCTION 


N-type materials have electrons as the majority carriers while p-type materials 
have holes as the majority carriers. What is important to semiconductor physics is what 
happens when these two types of materials come together [Ref. 2]. As can be seen from 
Figure | below, when the two materials come in contact, some electrons move across the 
junction to fill holes on the p-side and some holes move across the junction to the n-side 


[Ref. 2]. 


Depletion 
Region 





P-type N-type 
Figure 1. —p-n Junction 


When these electrons and holes migrate across the junction, an electric field 
forms. This area is referred to as the space charge region or depletion region, and it cre- 
ates a built in potential in the material. The energy band diagram in Figure 2 [Ref. 2] 
helps to give a better idea what happens when n-type and p-type materials are brought 


together. 
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Figure 2. Energy Band Model of a p-n Junction [From Ref. 2.] 


The Fermi level is constant throughout the material, so a straight line is drawn and 
the two band diagrams are shifted to match it. In order to have continuous conduction 
and valence bands in the material, the bands of the p- and mtype material must bend to 
meet each other [Ref. 2]. This bend in energy bands creates a built-in potential in the ma- 
terial equal to V,, as seen in Figure 3 [Ref. 2]. It is this potential that causes the voltage 


drop across a diode [Ref. 2]. 





Figure 3. —_ Built-In Potential [After Ref. 2.] 


Figure 4 shows how the junction reacts when a bias is placed across it [Ref. 2]. 
Part (a) shows the band diagram under equilibrium conditions with no bias present. This 


diagram was previously explained. 

Part (b) shows the junction under forward bias conditions. This is when the anode 
is held at a higher potential than the cathode. As more voltage is applied across the junc- 
tion, it acts against the built in potential V,,. Once the bias across the junction reaches a 
level that is equal to V,,, current will begin to flow through the junction. The diode will 


then act like a short circuit with a constant voltage drop of V,,. 


Part (c) shows the band diagram of the p-n junction under reverse bias conditions. 
This occurs when the voltage at the cathode is higher than the voltage at the anode. In 
this case, the potential across the junction acts to increase the level between the energy 


bands of the two materials thereby preventing current from flowing. 
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(c) Reverse Bias (V,<0) 
Figure 4. p-n Junction Under Bias [After Ref. 2.] 


C. PHOTOVOLTAIC CELL 


A photovoltaic cell [Ref. 3], shown in Figure 5, produces power using photogen- 
eration. When a photon of light with energy higher than the bandgap energy strikes the 
material, an electron-hole pair is created. The cell itself is simply a p-n junction which 
was discussed in the previous section. The electric field in the depletion region at the 
junction attracts the free electrons towards the cathode, the top layer in the figure, and 


attracts the holes towards the anode, the bottom layer. 


Light Generales Ligm Is Absorbed 
Electron and Hole at Back 
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Metal Contact 


© Electron 
© Hole 


Figure 5. Simple Solar Cell [After Ref. 3.] 


Current is carried away from the cell by contacts which are connected both on the 
top and bottom of the cell. The contact on the top of the cell is typically made up of a 
thin metal electrical grid to avoid shading the cell. The bottom metal contact, on the 
other hand, spans the entire surface of the cell. These contacts are connected to the ex- 


ternal load. 


D. INCIDENT SPECTRUM 


Solar cells can be designed to be optimized for different spectra based on their in- 
tended use [Ref. 3]. Cells intended for outer space are designed to the spectrum of Air 
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Mass Zero (AMO). This is the spectrum from the sun measured just outside the earth’s 
atmosphere. Figure 6 shows the AMO spectrum in | w/ m* / pm | . The data was meas- 
ured by The National Renewable Energy Laboratory (NREL) [Ref. 4]. The total power 


in this spectrum is approximately 1367 w/ m. 
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Figure 6. AMO Spectrum vs. Wavelength [Data for Plot is From Ref. 4.] 


As mentioned previously, the energy of a photon required to cause photogenera- 
tion in a material must be greater than the bandgap energy of the material. Because of 
this relationship, it is often useful to see the AMO spectrum as a function of energy in- 


stead of wavelength. The relationship of photon energy to wavelength is expressed as 


E=—, (2.1) 


(2.2) 


where h is Planck’s constant, c is the speed of light, and A is the wavelength of light 
[Ref. 3]. Figure 7 shows the AMO spectrum [Ref. 4] as a function of energy. 
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Figure 7. | AMO Spectrum vs. Photon Energy [Data for Plot is From Ref. 4.] 


As the sun’s spectrum passes through the earth’s atmosphere some wavelengths 
are absorbed and others are reflected. This causes the radiance at the earth’s surface to be 
a subset of the AMO spectrum. Other common spectra are Air Mass One (AM1) and Air 
Mass One Point Five (AM1.5). AM1 is the spectrum from the sun at earth’s surface on 
the equator when the sky is clear. Moving farther from the equator to higher latitudes 
and the presence of clouds and other weather conditions further reduces the radiance of 


the spectrum that reaches earth’s surface. AM1.5 is commonly used for testing terrestrial 
based solar cells [Ref. 3]. 
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E. POWER OUTPUT 


As discussed previously, a photovoltaic cell is simply a p-njunction. It therefore 
has associated with it the same current voltage relationship as a diode when there is no 


light incident on the cell. 


Figure 8 shows currents that contribute to the total current from a solar cell [Ref. 
5]. The curve in part (a) shows the photocurrent from the cell. This is a constant current 
level that is dependent on the amount power incident on the cell. The dark current curve 
in part (b) is an exponential curve that is a characteristic of the material. To find the cur- 
rent-voltage relationship of the solar cell, the dark current must be subtracted from the 


photocurrent. Part (c) of Figure 8 shows the resulting /-V curve of the cell. 


v Vv v 


(a) Photocurrent (b) Dark Current (c) |-V Curve 


Figure 8. Current in a Solar Cell [After Ref. 5] 


Figure 9 shows an example of a solar cell -V curve. The points on this curve of 
most importance are labeled. When the anode and cathode are shorted, the maximum 


current, called short circuit current, flows through the circuit [Ref. 3]. This is labeled as 
T,, on the figure. If the circuit is opened, no current will flow; however the maximum 


voltage, called open-circuit voltage, can be measured across the cell [Ref. 3]. 
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Figure 9. =‘ /-V Curve 


While the cell can operate anywhere on this curve, it is most efficient when oper- 
ated at the maximum power point, labeled P_,, in the figure [Ref. 3]. This is the point on 


the curve where the product of the operating voltage and current is the highest. The effi- 


ciency of the cell is calculated as [Ref. 3] 
n=—=, (2.3) 


where P, is the power in the incident spectrum; for AMO this is 1367 Ww/ m° [Ref. 3]. 


Another important characteristic of a solar cell is the fill factor, which indicates the 
squareness of the /-V curve [Ref. 3]. This is calculated as [Ref. 3] 
V, To. 
FrF=— | (2.4) 
Voi 


oc SC 


The dotted box that intersects 7,. and V,. in the figure represents the optimal /-V 


curve with a fill factor of one. The fill factor is a measure of how much of this dotted 


box the actual /-V curve occupies [Ref. 3]. 
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F. EFFICIENCY IMPROVEMENTS 


There are many factors that can be changed to help improve the efficiency of solar 


cells. 


1. Reflection 

An untreated solar cell can reflect as much as 36% of the incident light [Ref. 3]. 
Antireflective coatings such as silicon nitride can reduce the reflectivity of the cell to as 
low 5% [Ref. 3]. In addition to adding antireflective coatings, texturing the cell surface 


also reduces the amount of light reflected at the surface of the cell [Ref. 3]. 


Figure 10 shows an example of surface texturing [Ref. 6]. This cell was able to 
operate at 24% efficiency which is very high for a silicon cell [Ref. 6]. The inverted 
pyramids labeled in the figure are at such an angle so that light that is reflected from the 
surface of the cell, is reflected in such a way that it will hit one of the other pyramids 
where it may be absorbed [Ref. 6]. The high absorption of this cell led to its record set- 


ting efficiency. 





rear contact oxide 


Figure 10. Anti-Reflective Coating and Textured Cell Surface [From Ref. 6.] 
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Zi Low Energy Photons 

Only photons with energy higher than the bandgap energy of the material can cre- 
ate an electror hole pair in the cell. Photons with too little energy to create an electron 
hole pair will be absorbed and increase the temperature of the cell [Ref. 3]. This in turn 
reduces the efficiency of the cell. Some designs utilize optical filters to block low energy 


photons to avoid this problem [Ref. 3]. 


>. High Energy Photons 

Photons with energy significantly higher than the bandgap energy can still create 
an electrormhole pair, but the excess energy is given off as heat in the cell, thus reducing 
the efficiency [Ref. 3]. Some of these photons may simply pass through the cell without 
creating an electron-hole pair. A back surface reflector can reflect these high energy pho- 
tons back through the cell where they once again have the opportunity to create an elec- 


tron-hole pair [Ref. 3]. 


4. Cell Resistance 

There are several sources of resistance in a solar cell. There is a resistance in the 
material itself that impedes the path of free electrons and holes to the contacts [Ref. 3]. 
There are also junction resistances where the electrical contacts are joined to the cell 
[Ref. 3]. Finally there is resistance in the contacts [Ref. 3]. All of these resistances cause 


losses which can be significant especially in high current applications. 


5. Material Defects 

Defects in the crystal structure aid in trap assisted recombination which reduces 
cell efficiency [Ref. 2]. This is a significant source of power loss in inexpensive poly- 
crystalline and amorphous cells reducing efficiency by as much as 15% [Ref. 3]. There 
are also a significant number of material defects when lattice mismatch is a factor such as 
in multijunction cells. The effects of this in single junction cells can be reduced by using 
high quality crystals [Ref. 3]. The fewer the defects in the crystal structure, the higher 
the efficiency of the cell. In multijunction cells, a great deal of research has been done 


with ternary and quaternary compounds (compounds with three or four elements respec- 
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tively) [Ref. 3]. By altering the molar concentrations of the elements in these com- 
pounds, the lattice constant can be set to a desired value reducing the amount of lattice 


mismatch. 


6. Shading 

In order to collect free carriers from the cell, there must be metal contacts on the 
surface. These contacts, while providing a necessary function, also shade the parts of the 
cell they cover from incident light [Ref. 3]. Larger and more closely spaced contacts 
have the benefit of being able to collect more generated carriers and having lower internal 
resistance [Ref. 3]. Conversely they shade more of the cell causing fewer carriers to be 
generated. Smaller and more loosely spaced carriers shade less of the cell, but also have 
higher resistance [Ref. 3]. In addition when contacts are more loosely spaced on the sur- 
face of the cell, carriers are more likely to recombine before they can be collected by the 
contacts [Ref. 3]. By trial and error, it has been found that the optimal contact arrange- 
ment consists of equally spaced contacts covering approximately 8% of the cell’s surface 


[Ref. 3]. 


This chapter reviewed the basic principles of solar cells. The input spectrum and 
output power were covered along with methods for improving the efficiency. The next 


chapter will introduce a related concept, the thermophotovoltaic cell. 
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IW. THERMOPHOTOVOLTAIC CELLS 


This chapter discusses the basic principles of thermophotovoltaic (TPV) cells. 
Thermophotovoltaic cells are designed to convert the infrared (IR) region of the electro- 
magnetic spectrum into electricity. The goal of this thesis is to model and optimize these 


cells, so it is important to know some background information on them. 


A. HISTORY 


Although the ‘Father’ of thermophotovoltaics is unknown, the field can be traced 
at least as far back as 1956 to Dr. Henry H. Kolm [Ref. 7]. While working at MIT’s Lin- 
coln Laboratory, Dr. Kolm built a primitive TPV device using a Coleman camping la n- 
tern as a radiator and silicon solar cells for photoconverters. Although his system only 
produced one watt of power, it was a start. Throughout the 1960’s, researchers at MIT 
continued research on TPV devices focusing on germanium cells [Ref. 7.]. In the 1970’s 
thermophotovo lItaics was studied as a way to provide a portable power source for the U.S. 
Army [Ref. 7]. Early TPV generators consisted of standard silicon or germanium solar 
cells surrounding a heat source [Ref. 7]. Around 1980 a general interest in the area was 
lost because of the exceedingly low efficiency of these early TPV cells. With the recent 
advancements in III-V semiconductor technology (materials made of a combination of 
group III and group V elements from the periodic table), there was a possibility for sig- 
nificantly higher efficiency cells [Ref. 7]. Recent research has shown that gallium anti- 
monide (GaSb) and indium gallium arsenide (InGaAs) are two of the most promising ma- 


terials for new TPV cells. 
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B. INCIDENT SPECTRUM 


1. Blackbody Radiation 
Objects with a temperature above zero Kelvin radiate energy. In order to estimate 
the amount of energy radiated, it is convenient to introduce the concept of a blackbody. 


A blackbody is characterized by three key properties [Ref. 8]: 
1. 100% absorption of all incident radiation 
2. Maximum possible emission at any temperature and wavelength 


3. Diffuse emitter, radiates energy equally in all directions 


An ideal emitter, a blackbody, radiates energy over the whole electromagnetic 
spectrum as a function only of temperature and wavelength. A blackbody emitter radi- 
ates energy per wavelength into a unit solid angle in a distribution described by Planck’s 


Law for spectral radiance as [Ref. 8] 


(3.1) 


La) WwW |- 2hc? 


m’ -sr-um 


where /: is Planck’s constant, c is the speed of light in a vacuum, and k is Boltzmann’s 
constant. T and A are the only two variables and refer to the blackbody temperature in 
Kelvin and the wavelength of energy emitted in meters. The units of spectral radiance 


are bracketed. 


As can be seen from Planck’s Law for spectral radiance, at a given temperature, 
the amount of radiation varies as function of wavelength. An important feature is that at 
any wavelength, the magnitude of radiant energy increases as temperature increases [Ref. 


8]. This means that a hotter blackbody will emit more energy than a colder one. 


Another important feature of Planck’s Law is the distribution of radiant energy. 
More radiant energy is emitted at shorter wavelengths for hotter blackbodies [Ref. 8]. 


This phenomenon is described by Wien’s Displacement Law [Ref. 8], 
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This equation gives the wavelength at which the maximum radiation is emitted 
from a blackbody at a given temperature. It can easily be seen that the location of this 


peak is inversely proportional to the temperature of the blackbody. 


In order to find the total power emitted by a blackbody, it is necessary to integrate 
Planck’s Law over all wavelengths. This calculation yields the Stefan Boltzmann Law 
[Ref. 8], 

EAGT, (3.3) 
where 6 =5.670x10* [ w/ m’/ K*] . Here it can be seen that the total radiated power is 


dependent only on temperature. 


When working with solar cells, the spectrum of interest is AMO. This spectrum 
corresponds to the radiation from blackbody at a temperature of 5800K. Because energy 
is radiated outward from a body at an average angle perpendicular to its surface, the mag- 
nitude of radiation in a given direction is equal to only one-fourth the total radiation [Ref. 
8]. This corresponds to 7 steradians. With this information, the equation for the hemi- 
spherical spectral radiant flux in any direction from a flat surface emitter can be derived 
as [Ref. 8] 
2m he? 


x [Fl] 


In order to find the AMO spectrum, there is one more factor that must be consid- 


un] 


m?- um 


=nL(T)= (3.4) 


ered. The magnitude of radiation emitted from a source is inversely proportional to the 
square of the distance from the source. The factor used to calculate the AMO spectrum is 
found by taking the square of the ratio of the sun’s radius to the mean distance between 


the earth and the sun. This factor, referred to as the dilution factor, comes out to 


2.16510” [Ref. 8]. The final equation for the AMO spectrum can now be expressed as 


WwW f (2nhc’) 
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where /f is the dilution factor. 


Using this equation, Figure 11 below shows the spectrum from a 5800K black- 
body adjusted by the dilution factor. It is plotted on top of the actual AMO spectrum as 
measured by NREL [Ref. 4]. It can be seen that the two spectra are very closely corre- 


lated in shape, magnitude, and area. 
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Figure 11. AMO vs. 5800K Blackbody Spectrum 


From Planck’s distribution law, the magnitude of radiation at a given wavelength 
increases with temperature, so the magnitude of radiation from the sun is much larger 
than the magnitude of radiation from any heat source with a lower temperature. Unlike 
photovoltaics, however, where the radiant source, the sun, is millions of miles away, a 
heat source used for thermophotovoltaic applications can be placed in very close prox- 
imity to the cell. Although the sun radiates much more energy than a heat source, be- 
cause of the distance between the source and the cell, the intensity of radiation at the cell 


is much higher for a heat source. 


20 


Zz Radioisotope Thermoelectric Generator Spectrum 

When dealing with solar cells, the input spectrum is a subset of the spectrum from 
the sun [Ref. 3]. With TPV cells, the input spectrum is dependent on the temperature of 
the radiator. Table 1 shows several temperatures and the wavelength at which their radia- 


tion spectrum is a maximum as calculated from Equation 3.2. 


T(K) | ? ca (umn) 
300 | 9.66 

10002507 
1300 [223 
2000; 1.45 

SE 
5000 | 0.58 

5800) 0.50 | 
7000; O41 
100007 0.39 


Table 1. Temperature vs. Peak Wavelength 


The previous section discussed blackbody radiation and compared a 5800K black- 
body radiation curve to the AMO spectrum in Figure 12. Because these two were so 
closely correlated, it was assumed that a blackbody curve would be a good best case ap- 
proximation for the input spectrum from a radioisotope thermoelectric generator (RTG) 


heat source. 


Current RTG’s use Pu-238 as a heat source. In a pure form, Pu-238 can reach 
surface temperatures of 1300K [Ref. 9]. For simulations in this research, the input spec- 
trum was assumed to be a 1300K blackbody. Figure 24 shows the spectrum that was 


used in simulations. 
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1300 K Blackbody Spectrum 
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Figure 12. 1300K Blackbody Spectrum 


Cc. THERMOPHOTOVOLTAIC SYSTEM 


There are more parts to a thermophotovoltaic system than just the TPV cell. Fig- 


ure 13 shows the four main components of a typical TPV system. 
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Figure 13. TPV System 
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1; Heat Source 
The first component is a heat source. The temperature of the heat source is a very 
important factor in the efficiency of the system. This relationship will be discussed in 


more depth in the next section. 


2. Selective Emitter 

The second component in the system is a selective emitter. The purpose of this 
emitter is to help shape the radiation spectrum from the heat source. As with solar cells, 
radiation absorbed by the cell that does not have enough energy to create an electron hole 
pair just adds heat to the cell, thereby lowering the efficiency. Rare earth oxides such as 
ytterbium oxide (Yb2O3) and yttrium aluminum garnet (YAG) have proven to be very 
effective materials because they emit radiation in a line spectrum instead of a continuous 
spectrum like a heat source but they are currently in the experimental stages and are not 
always found in TPV systems [Ref. 10]. A selective emitter was not used in these simu- 
lations because the response of the materials used was desired over a wide range of the 
spectrum. This data is important when designing multijunction cells which will be dis- 


cussed later. 


3. Filter 

The filter performs a similar function to the selective emitter. It is designed to re- 
flect incident radiation that does not have enough energy to create an electror-hole pair in 
the TPV cell. The radiation is reflected back to the emitter. This helps to keep the tem- 
perature of emitter as high as possible which improves the efficiency of the system [Ref. 
11). 


Figure 14 shows an experimental resonant gold film filter [Ref. 11]. Designed by 
Energetics Development Technology (EDTEK), this filter is etched from a solid gold 
film in a pattern designed to resonate at a desired frequency [Ref. 11]. The filter can be 
tuned to resonate at any desired frequency depending on the intended spectrum. Tests 


have shown that the filter has very good transmission in the desired band and reflects 
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over 98% of out-of-band radiation [Ref. 11]. Filters were not included in the models 


built here for the same reasons that selective emitters were not used. 
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Figure 14. | Resonant Filter [After Ref. 1 Ld 


Figure 15 shows the frequency response of a filter built with this technique along 
with the AMO spectrum [Ref. 11]. This particular filter was designed to resonate at two 
distinct frequencies. This capability is useful for multijunction cells which consist of in- 


dividual cells that are designed to respond to different wavelengths of light. 
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Figure 15. Resonant Filter Frequency Response [From Ref. 11.] 
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4. Photo-Electric Converter 

Conventional RTG’s use a thermocouple to convert heat into electricity [Ref. 13]. 
A thermocouple produces power using what is known as the Seebeck Effect. Seebeck 
discovered that when two materials from the thermoelectric series were connected, a po- 
tential could be measured across them if the two materials were at different temperatures 
[Ref. 12]. In addition, if the two materials were connected with a wire so as to close the 
circuit, a small current was forced to flow. In an RTG, a SiGe thermocouple is placed in 
the vicinity of the heat source with one side exposed to the heat and the other side to a 
heat sink [Ref. 13]. The electrical power produced is proportional to the difference in 
temperature between the hot and cold junctions. These devices do a relatively poor job of 
power conversion, however, only achieving efficiencies of a few percent in optimum 


working conditions [Ref. 12]. 


This thesis will explore the thermophotovoltaic cell as a candidate for replacing 
the thermocouple in RTG’s. Such a cell operates in a manner very similar to how a solar 
cell works. The difference between the two is the band of the electromagnetic spectrum 
the cell is designed to respond to. Where a solar cell responds to the visible band of the 


spectrum, thermophotovoltaic cells respond to the infrared region. 


In order for the TPV cell to respond to radiant energy from an RTG source, the 
energy of the radiant photons must be greater than the energy bandgap of the material. 
As previously mentioned, photon energy is inversely proportional to the temperature of 
the radiant body. From Table 1 it can be seen that the radiation from a 1300K blackbody 
source peaks at 2.23 um which is a much longer wavelength than the 0.45 um the sun’s 
spectrum peaks at. This means that materials used for TPV cells must have very low en- 
ergy bandgaps in order to respond to such long wavelengths of light. Table 2 shows sev- 
eral different semiconductor materials with their respective energy bandgaps and the 


maximum wavelength of light they can respond to as calculated from Equation 2.2. 
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GaAs ; 1.42 ; 0.87 
Ss 
Gass 072 1.72 
GSO 066 188 
“inGasb | 0.6 12.07 
InGaAsSb | 0.55 | 2.23 
InGaAs | 0.55 | 2.23 
TnAssbP 0.30 | 3.18 


Table 2. |. Maximum Frequency Response of Various Materials 


In order for the peak of the 1300K blackbody spectrum to occur in the energy 
band that a TPV cell can respond to, the cell must be made of a material with an energy 
bandgap lower than 0.55 eV. As can be seen from Table 2 there are materials that have 
such a low bandgap; however there are tradeoffs to be made. The opercircuit voltage of 
a cell is closely related to the energy bandgap of the material. As the energy bandgap 
goes down, so does the oper-circuit voltage. A material with a bandgap of 0.39 eV may 
only produce an open-circuit voltage of just a few tens of millivolts or less. Such a low 
voltage drastically reduces the maximum power from the cell in turn reducing the effi- 


ciency. 


Gallium antimonide (GaSb) is the most commonly published material used for 
TPV cells. Although it does not have a low enough energy bandgap to respond to most 
of the radiant energy from an RTG heat source, the cells are relatively easy to fabricate 
and are capable of producing a significant amount of power per square centimeter. In- 
dium gallium arsenide (InGaAs) is another popular material for TPV cells. As a ternary 
compound, many of its electrical and optical properties including energy band gap 
change with varying molar concentrations of indium and gallium. Asa 0.55 eV cell, its 
low operrcircuit voltage is not very suitable for single junction cells, but it is attractive 
for use in multijunction cells [Ref. 14]. For these reasons, GaSb and InGaAs are the ma- 


terials that will be studied in this thesis. 
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This chapter gave an in depth discussion of thermophotovoltaics including the 
cells themselves and the spectrum they are designed to respond. The next chapter will 
cover the systems that thermophotovoltaic cells are designed to operate in, radioisotope 


thermoelectric generators. 
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IV. RADIOISIOTOPE THERMOELECTRIC GENERATOR 


This chapter presents the radioisotope thermoelectric generator (RTG). Asa 
source of space nuclear power, the use of RTG’s, while proven safe, reliable, and neces- 
sary for deep space missions, has been greatly debated. The history of RTG’s will be 
discussed including an account of their successes. Environmental issues regarding the 
use of RTG’s will also be mentioned along with the many measures in place to maintain 
their impeccable safety record. Finally this chapter also includes about possible im- 


provements to RTG’s necessary for higher efficiency operation. 


A. HISTORY 


As space missions have begun to take man-made vehicles and satellites farther 
and farther from the sun, the need for a better means of power generation has become 
more important. The intensity of light from the sun reduces with the square of the dis- 
tance from it meaning that massive solar arrays would be needed in order to produce 
minimal power on distant missions. NASA has calculated that on missions such as Cas- 
sini, the weight of solar arrays needed to provide the desired amount of power would 


more than double the weight of the entire satellite [Ref. 3]. 


An alternative to these massive arrays is to bring a fuel supply that will be capable 
of providing adequate power over the entire lifespan of the vehicle. It is not feasible to 
use just batteries because the number of batteries required for such a long duration would 
be far too heavy [Ref. 3]. Conventional fuels are also not viable options due to the poor 


mass to power ratios [Ref. 3]. 


Since 1961, the United States Department of Energy (DOE) has supplied NASA 
with RTG’s powered by Pu-238 [Ref. 13]. Table 3 lists some missions that have made 
use of these RTG’s [Ref. 15]. 
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Misson Type 


Status 


Operated for 15 years. 
Shut down but still 


operational. 


RTG operated as 
planed. Electrical 
problems caused 


failure after 9 months. 


RTG operated for 6 
years. Navigation failed 


after 1.5 years. 


Mission aborted. RTG 
burned up on reentry 


as designed. 


Successfully achieved 


orbital operations. 


Mission aborted. 
RTG recovered and 


recycled. 


RTG's operated for 


more than 2.5 years. 


RTG shut down after 


8 years of operation. 


Mission aborted. Heat 


Source returned. 


RTG shut down after 
6.5 years of operation. 
~RTG shut down after _ 
6 years of operation. 
~~ RTG still operation. — 

Spacecraft now 
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Spacecraft 


Apollo 16 


Apollo 17 


Pioneer 11 


Viking 2 


Voyager 2 


Voyager 1 


Galileo 


Ulysses 


Cassini 
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Lunar Surface 


Planetary 
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Status 


RTG shut down after 


5.5 years of operation. 


RTG shut down after 


5 years of operation. 


Spacecraft operated to 


Jupiter and Saturn. 


RTG shut down after 
6 years of operation. 


RTG operated 4 years 


until relay link was lost. 


TG still operating. 


ve) 


Operated to Jupiter, 


N 


aturn, Uranus, Neptune. 


“RIG still operating. — 
Operated to Jupiter, 
Saturn, and beyond. 

~RTG still operating. — 
Spacecraft orbiting 
around Jupiter. 

“RIG still operating. — 
Spacecraft in polar 
orbit around the sun. 


RTGs still operatin. 


Spacecraft en route 
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to Saturn. 


. _RTG Missions [After Ref. 15.] 


This table shows that RTG’s are capable of providing power over long periods of 


time. The first RTG sent into space in June of 1961 on the Transit 4A spacecraft is still 
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providing power [Ref. 15]. Figure 16 shows a rendering of a current RTG module [Ref. 
16]. 
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Figure 16. .RTG Module [From Ref. 16.] 


B. FUEL SUPPLY 


An RTG fuel supply must be able to balance a high power-to- mass ratio with a 
long lifetime and reasonable safety. As mentioned previously, current RTGs have used 
Pu-238 as a heat source. Pu238 has a half-life of 87 years, meaning that after 20 years of 
operation it will still produce 80% of the power it produced when new [Ref. 17]. Figure 
17 shows the power output versus time for a Pu:238 heat source [Ref. 17]. This isotope, 


when pure, has a relatively large power to mass ratio of 540 W/kg [Ref. 17]. 
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Figure 17. Pu-238 Power Output [From Ref. 17.] 


Other fuel supplies have been considered for use with RTG’s as well. Strontium 
90 (Sr-90) for example was considered as a possible fuel source. With a half life of 28.1 
years, Sr-90 decays much faster than Pu-238 [Ref. 18]. After 20 years, it will only pro- 
duce approximately 60% of the power it originally did [Ref. 18]. In addition, the power 
to mass ratio when new is only 460 W/kg which is more than 15% lower that that of Pu 


238 [Ref. 18]. 


Although the use of Pu:238 has startled many people, it is actually a rather safe 
fuel. In order to demonstrate to the public how safe the use of Pu-238 is, President Et 
senhower had the picture in Figure 18 published in the newspaper taken in the Oval Of- 
fice at the White House with senior members United States Atomic Energy Commission 


(AEC) and an early RTG module [Ref. 17]. 
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Figure 18. “President Shows Atom Generator” [From Ref. 17.] 


C. SAFETY 


if Accident History 

In total the United States Department of Energy has provided 44 RTG’s and more 
than 240 general-purpose heat source (GPHS) units since 1961 with only three accidents 
[Ref. 17]. The first of these accidents occurred on the Transit 5-BN-3 launched in April 
of 1964. This mission had to be aborted when the launch vehicle failed. The vehicle and 
RTG burned up in the atmosphere upon reentry where the Pu-238 fuel source was vapor- 
ized and dispersed into the upper atmosphere as it was designed to do [Ref. 17]. Shortly 
after this accident, the RTG module was redesigned to survive atmospheric reentry with 


the fuel source intact [Ref .17]. 


The second accident occurred on the Nimbus-B-1 launched in May of 1968 [Ref. 
17]. This mission had to be aborted shortly after launch by a range safety destruct [Ref. 
17]. The newly designed RTG module held up through reentry and was recovered with- 


out any radiation leakage. The fuel source was reused on later missions [Ref. 17]. 


The third and final accident occurred on Apollo 13 launched in April of 1970. 
This mission carried an RTG intended to provide power to a lunar seismic station [Ref. 
17]. The RTG was located on the lunar module which broke up on reentry into earth’s 
atmosphere [Ref. 17]. The RTG remained intact and is now located at the bottom of the 
Tonga Trench in the Pacific Ocean [Ref. 17]. 
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2. Safety Measures 
Many features have been added to RTG’s to ensure their safe operation under all 


conditions. Figure 19 shows a diagram of a modern RTG unit [Ref. 19]. 
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Figure 19. RTG Safety Measures [From Ref. 19.] 


It can be seen from the figure that the fuel pellet makes up a very small percent- 
age of the unit while added safety features occupy the majority of the area. The fuel itself 
has even been modified for safety purposes. The original Pue238 metal has been replaced 
with plutonium dioxide [Ref. 17]. This heat resistant ceramic form prevents the fuel from 
vaporizing under extreme heat conditions such atmospheric reentry [Ref. 17]. The mate- 
rial is also insoluble and has a very low chemical reactivity making it safe even if it ends 
up in the ocean like the module from Apollo 13 [Ref. 17]. This ceramic form also influ- 
ences how the material breaks up. Instead of becoming a fine powder that can easily 
spread, it fractures into large nomrespirable chunks [Ref. 17]. As an added safety meas- 
ure, the fuel is separated into small pellets that are each individually encapsulated in a 
heat shield and impact shell. The module itself is enclosed in an iridium capsule and sur- 
rounded by graphite blocks [Ref. 17]. These are high strength and corrosion and heat re- 
sistant materials. The RTG unit is designed to withstand conditions it could never realis- 


tically encounter. 
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The Pu-238 in RTG’s is an alpha emitter [Ref. 17]. The majority of alpha parti- 
cles can be absorbed by a piece of paper or the epidermal layer of human skin [Ref. 17]. 
RTG’s do not utilize either fission of fusion in producing heat [Ref. 17]. Because of this, 
an RTG could never explode like a nuclear bomb. In addition, no accident involving an 
RTG could ever create the acute radiation sickness associated with nuclear explosions or 


meltdowns [Ref. 17]. 


a Safety Testing 
As previously mentioned, RTG modules are designed to endure conditions they 
could never realistically face. In order to ensure the safety of the modules, they are tested 


under extremely harsh conditions [Ref. 17]. 


a. Fire 
Extreme heat upon atmospheric reentry and other unforeseen accidents 
can cause fires onboard a spacecraft. RTG’s have been directly exposed to intense fires 


with minimal damage [Ref. 17]. Most importantly, no fuel was released during exposure. 


b. Blast 

A shuttle explosion can cause tremendous shockwaves. The RTG module 
was tested under blast waves stronger than those expected from such an explosion and 
once again remained intact without releasing any nuclear fuel or harmful radiation from 


the unit [Ref. 17]. 


c. Reentry 

In the event that a satellite degrades in orbit enough to reenter the earths 
atmosphere, it will experience intense amounts of heat. An RTG is designed to survive 
reentry while preventing the Pu-238 fuel source from vaporizing in the atmosphere. This 
environment was tested using an arc-jet furnace, and the RTG was able withstand the 


heat safely [Ref. 17]. 
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d. Earth Impact 

Since RTG’s were redesigned to survive atmospheric reentry intact, they 
also had to be strong enough to survive impact with earth’s surface. RTG’s were tested 
at the maximum estimated velocity of the module approaching earth’s surface of 120 
miles per hour. These tests showed that in an impact with water, soil, and sand the mod- 
ule would remain safely intact. An impact with rock or concrete, however, sometimes 
caused slight leakage resulting in small amounts of low-level radiation in localized areas 


[Ref. 17]. 


e. Immersion in Water 

Because the vast majority of earth’s surface is covered in water, RTG’s 
must be able to survive in an aquatic atmosphere. Long term exposure to sea water 
proved the protective iridium capsule to be corrosion resistant and the fuel to be highly 


insoluble [Ref. 17]. 


f Shrapnel 

In the event of an explosion or impact with space debris, an RTG could be 
subjected to shrapnel. Aluminum and titanium bullets were used to simulate such shrap- 
nel from a launch explosion. Fragments at speeds much greater than those predicted to 


be present in such an event were unable to penetrate the module [Ref. 17]. 


g. Large Fragments 

In the event of a solid rocket booster failure, the RTG module could be 
subject to impact by large fragments of material. To simulate this scenario, large steel 
plates were fired at an RTG. The studies showed that damage only occurred in the 
unlikely event that the corner or sharp edge of a large fragment struck the module [Ref. 


17]. 
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D. FUTURE IMPROVEMENTS 


Current RTG’s use silicon/germanium thermocouples discussed in the previous 
chapter to convert thermal power from the Pu-238 heat source into electrical power [Ref. 
13]. These thermocouples operate at very low efficiency levels of only a few percent. 
The key to improving RTG power production is to find a more efficient way to convert 


thermal power into electrical power. 


1. Dynamic Isotope Power Systems 

Dynamic isotope power systems (DIPS) convert thermal energy into electrical en- 
ergy using moving parts. Created in 1816 by Robert Stirling, the Stirling engine is the 
most important dynamic isotope power system to space power systems [Ref. 20]. Unlike 
combustion engines, Stirling engines have a closed gas chamber [Ref. 20]. There are no 
explosions that take place in a Stirling engine and there is no gas exhausted from the en- 
gine. The power to run this engine comes from an external heat source [Ref. 20]. Figure 


20 shows a simplified example of a Stirling engine [Ref. 20]. 


om ye 





Figure 20. Simplified Stirling Engine [From Ref. 20.] 


There are four steps in the Stirling cycle [Ref. 20]: 


1. Heat is added to the hot cylinder (left), causing pressure to build. This pushes 


the left piston down. 
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2. The left piston moves up and the right piston moves down. This pushes the 
heated gas from the left cylinder into the cold cylinder on the right where the gas is 
quickly cooled down. When the gas cools, the pressure goes down making it easier to 


compress in step 3. 


3. The piston in the cold cylinder on the right compresses the gas. Heat gener- 


ated by this compression is almost instantly removed by the cooled cylinder. 


4. The right piston moves up while the left piston moves down forcing the gas 
into the heated cylinder on the right. This cylinder quickly heats the gas raising the pres- 


sure in the cylinder. 


The cycle then repeats. 


In this four-step cycle, power is only created in the first step. There are two ways 
to increase the power output of a Stirling engine. The first method is to increase the 
power output of stage one [Ref. 20]. This is accomplished by raising the pressure of the 
gas in the hot cylinder which is done by increasing the temperature of the external heat 
source. The second method for increasing the power output is lowering the gas pressure 
in the cold cylinder [Ref. 20]. This is accomplished by lowering the temperature of the 
cylinder. The Stirling engine, like a thermocouple, produces a power output related to 


the difference in temperature of the hot and cold cylinders. 


The Stirling engine is not suitable for widespread use for two main reasons. Be- 
cause the heat source is external to the engine, it takes some time for the engine to warm 
up before it can produce useful power [Ref. 20]. Secondly, the engine is slow to respond 


to changes in the power supply [Ref. 20]. 


Stirling engines do however, show promise as a means of power generation for 
space based systems. NASA Glenn has been developing Stirling engine technology since 
the 1970’s and has been developing radioisotope power systems for deep space missions 
since 1990 [Ref. 21]. Figure 21 shows a Stirling generator developed by NASA Glenn 
[Ref. 21]. 
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Figure 21. 110 Watt Electric Stirling Generator [From Ref. 21.] 


This figure shows two synchronized 55 watt electric generators. This configura- 
tion reduces vibrations in the system. The addition of internal precision springs and 
dampers causes this engine to operate at a resonant frequency that increases the effi- 
ciency of the system. This model has been tested to work at an efficiency level of over 
20% [Ref. 21]. 


2. Alkaline Metal Thermal-to-Electric Conversion 
Figure 22 shows the schematic diagram of a simple alkaline metal thermal to- 


electric conversion (AMTEC) system [Ref. 22]. 


Hh ays 





A-Heat Source B- Liquid Sodium C- Heat Sink 
D- Sodium Vapor E- Beta-Alumina Solid Bectralyte 


Figure 22. . AMTEC Schematic [From Ref. 22.] 
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The key to this converter is the beta-alumina solid electrolyte (BASE). This con- 
ducts positive sodium ions much better than it conducts sodium atoms or electrons [Ref. 
22]. When heat is applied to one side, a sodium pressure difference occurs across the 
base and the positive ions move through the base and collect on the low pressure side 
leaving electrons to collect on the high pressure side [Ref. 22]. This creates a potential 
across the base that can force a current to flow through a load [Ref. 22]. Similar to the 
previously discussed thermoelectric converters, the power produced by an AMTEC sys- 
tem is dependent on a large temperature difference between the hot and cold ends, but 


currently the efficiency of such systems is very low[Ref. 22]. 


3. Thermophotovoltaic Cells 

Thermophotovoltaic (TPV) cells are currently the most promising thermoelectric 
converters for replacing the thermocouple in RTG’s. Although they are currently not ca- 
pable of the operating efficiencies of Stirling engines, TPV cells have the advantage of no 


moving parts [Ref. 3]. These cells are discussed in more depth in the previous chapter. 


The radioisotope thermoelectric generator was introduced in this chapter. The 
current system as well as possible future improvements were covered. The next chapter 


will begin the discussion on the simulation software used in this thesis. 
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V. SIMULATION SOFTWARE 


There is currently a great deal of research being done on building accurate models 
of solar and thermophotovoltaic cells. While some models can provide relatively accu- 
rate results in limited cases, there is no widely available software package that is capable 
of modeling all of the necessary parameters needed to provide models of cells of various 
materials under many different spectra. This thesis explores the use of Silvaco Virtual 


Wafer Fabrication software to accomplish this task. 


A. SILVACO 


Silvaco International is a software company that provides programs for modeling 
all areas of electronics including digital and analog circuits [Ref. 23]. The company has 
software ranging from simple Spice modeling to advanced integrated circuit layout and 
extraction tools. The Virtual Wafer Fabrication package is part of Silvaco’s TCAD suite. 
It was designed to accurately emulate physical wafer manufacturing [Ref. 23]. The pack- 
age allows users to easily optimize all aspects of IC design before actually fabricating 
them. ATLAS is the main program in this Virtual Wafer Fabrication package that will be 
used in this research. This program can handle one-, two-, and three-dimensional designs 


and contains the following comprehensive models [Ref. 23]: 
e DC, AC small-signal, and full time-dependency 
e Drift-diffusion transport models 
e Energy balance and hydrodynamic transport models 
e Lattice heating and heatsinks 
e Graded and abrupt heterojunctions 
e Optoelectronic interactions with general ray tracing 
e Amorphous and polycrystalline materials 


e General circuit environments 
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e Stimulated emission and radiation 

e Fermi-Dirac and Boltzmann statistics 

e Advanced mobility models 

e Heavy doping effects 

e Full acceptor and donor trap dynamics 

e Ohmic, Schottky, and insulating contacts 

e SRH, radiative, Auger, and surface recombination 
e Impact ionization (local and nor local) 

e Floating gates 

e Band-to-band and Fowler-Nordheim tunneling 
e Hot carrier injection 


e Thermionic emission currents 


ATLAS is a powerful program, but it relies on the aid of several other Silvaco 
programs to build an accurate model of a circuit. The program Deckbuild is the run time 
environment for this software package. A TPV cell can be defined graphically from the 
command line in this program in a method that will be discussed later. ATHENA, the 
process simulator, actually builds the devices defined by the commands entered in the 
Deckbuild interface [Ref. 23]. DevEdit is another program that allows a user to define 
cell graphically by defining mesh lines and nodes on a plane [Ref. 23]. This can also be 
done directly through Deckbuild as will be discussed later. In addition to run time out- 
put, ATLAS also writes *.log files which contain the output data solved for during simu- 
lations. The output gathered from these simulations was voltage and current levels and 
photogeneration rates. The included program TonyPlot can input these *.log files and 
present the information in a graphical form. Figure 23 shows a block diagram of how 


these programs interact to solve a particular TPV cell [Ref. 23]. 
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Figure 23. Simulator Block Diagram [From Ref. 23.] 


B. DEFINING A CELL 


In this thesis, solar and TPV cells were defined graphically in ATLAS through the 
Deckbuild interface. In order to fully define a cell, it was necessary to enter not just di- 
mensions, but also material properties. The following format was used in building each 


model. 
e Defining constants (cell dimensions, doping levels, shading, etc.) 
e Defining mesh boundaries using the previously defined constants 
e Building the X- mesh 
e Building the Y-mesh 
e Defining regions (anode, cathode, contacts, etc.) 
e Placing electrodes 
e Defining doping levels in each region 
e Defining material properties (electrical, optical, band, etc.) 
e Adding incident spectrum 
e Solving 
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iF Defining Constants 

ATLAS allows a user the option of building a model directly by inputting dimen- 
sions in numbers or indirectly by entering dimensions as variables names that have been 
previously defined. In this thesis, it was necessary to be able to change various dimen- 


sions quickly and easily so the second method of defining a cell was utilized. 


The first portion of the code was dedicated to defining constants that would later 


be used in calculations. For example 
set cathodethickness=10 
set cathodedoping=1.000000e+020 
set anodethickness=.37 
set anodedoping=9.8e19 
set bsfthickness=15.000000 
set bsfdoping=1.000000e+020 


defines constants that represent the thickness and doping level for the cathode, anode, and 
BSF layers of a cell. These variables can now be used later in the program. The benefit 
to defining dimensions in this manner is that by simply changing the variable definition, 
the whole model is altered to reflect this change making iterations quick and easy to exe- 


cute. 


Zi Defining Mesh Boundaries 

ATLAS solves equations at many different nodes in a defined area. These nodes 
are defined by meshes. The constants defined in the first part of the ATLAS code are 
used in defining mesh boundaries. The following commands are an example of how 


mesh boundaries are defined in a cell. 
set bsfL=0 
set bsfH=$bsfL-$bsfthickness 


set bsfdivs=$bsfthickness/$divs 
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set catL=$bsfH 

set catH=$catL-$cathodethickness 
set catdivs=$cathodethickness/$divs 
set anoL=$catH 

set anoH=$anoL-$anodethickness 
set anodivs=$anodethickness/$divs 


The dollar sign is an indication to ATLAS that the following string is a user de- 
fined constant. In ATLAS, moving up on the Y-axis is in the negative direction so these 
commands define a cell by starting at Y =0 and building upwards including divisions in 


each region. 


Si Building the X- and Y- Meshes 
ATLAS requires that each region boundary fall on a mesh. In order to ensure that 
this is the case, mesh lines are manually drawn at each of these boundaries using the fol- 


lowing commands. 
x.mesh loc=0 spac=$cellwdiv 
x.mesh loc=$contLL spac=$capwdiv 
x.mesh loc=$contRR spac=$capwdiv 
x.mesh loc=$cellwidth spac=$cellwdiv 
y-mesh loc=$contH spac=$contdivs 


y-mesh loc=$capH spac=$capdivs 





y.-mesh loc=$catH spac=$catdivs 
y-mesh loc=$anoH spac=$anodivs 
y-mesh loc=$bsfH spac=$bsfdivs 


y.mesh loc=$bsfL spac=$bsfdivs 
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Each of these command lines defines both a location and a spacing interval. 
ATLAS will draw mesh lines at each defined location as well as lines at the defined spac- 
ing interval. When two defined mesh lines with different spacing intervals are next to 
each other, ATLAS will automatically grade the mesh spacing between the two lines. 


Figure 24 shows an example of a typical mesh. 
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Figure 24. Typical Mesh 


When defining a mesh, it is important to take careful note of the spacing. In order 
to get an accurate solution, the mesh must be fine enough to get an accurate representa- 
tion of the cell. If the mesh is too fine, however, it will take an unnecessarily long time to 


solve all of the nodes. 


4. Defining Regions 

Each layer of a cell must be defined by a material and region boundaries. As pre- 
viously mentioned, these boundaries must fall on already defined mesh lines. The fol 
lowing commands are used to define regions in the GaSb cell shown in Figure 25. The 


bracketed numbers represent the value of each variable used in micrometers and are not 
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included in the actual code. This cell has a 0.1 um thick cap, 0.37 um anode, and a 
10 um thick cathode. The BSF layer defined in these commands cannot be seen in Fig- 


ure 25. 


#Anode 


region num=5 material=GaSb x.min=0 x.max=$cellwidth[S00] 


y-min=$anoH[-350.37] y.max=$anol[-350] 
#Cathode 


region num=6 material=GaSb x.min=0 x.max=$cellwidth[5S00] 


y-min=$catH[-350] y.max=$catL[-340] 
#BSF 


region num=7 material=GaSb x.min=0 x.max=$cellwidth[500] 


y-min=$bsfH[-340] y.max=$bsfL[0] 


Figure 25 shows a plot of cell regions generated by TonyPlot. 
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Figure 25. Cell Regions 
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5. Placing Electrodes 

Once the cell regions are fully defined, electrodes must be defined to electrically 
connect to the cell. ATLAS allows a user to define an ideal contact or specify a partic u- 
lar material with nor ideal properties. It is even possible to add external discrete compo- 
nents to the cell such as resistors, inductors, and capacitors. Electrodes are defined using 


the following commands. 


electrode name=anode material=Gold x.min=$contLL x.max=$contRR 


y-min=$contH y.max=$contL 


electrode name=cathode material=Gold x.min=0 x.max=$cellwidth 
y.min=$bsfL y.max=$bsfL 


6. Defining Doping Levels 
The doping level in each region can be defined independently using the following 


commands. 
# Anode 
doping uniform region=5 p.type conc=$anodedoping 
# Cathode 
doping uniform region=6 n.type conc=$cathodedoping 


This defines a constant doping level in a region. The boundaries between two re- 
gions will have abrupt junctions. For more realistic and more complex doping schemes, 


Silvaco provides a C-Interpreter, which allows a user to define any conceivable doping. 


7. Defining Material Properties 
Using the following commands, the user can define as many material properties 


as are available. 
material material=GaSb COPT=5e- 10 


material material=GaSb EG300=.726 PERMITTIVIT Y=14.4 
AFFINIT Y=4.06 


50 


material material=GaSb NC300=2.1e17 NV300=1.8e19 
material material=GaSb index.file=GaSb.opt 


These commands define optical recombination, energy band gap, permittivity, 
electron affinity, available conduction and valence band states, and optical parameters for 
the material gallium antimonide (GaSb). The index file “GaSb.opt’ contains refractive 
index and extinction parameters for a span of wavelengths in the area of interest. This is 
one of the most important material properties to accurately define and subsequently one 


of the most difficult to find. Figure 26 shows a plot of the data from ‘GaSb.opt’. 
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Figure 26. “GaSb.opt’ [Data for Plot is From Ref. 21] 


8. Incident Spectrum 


ATLAS allows a user to define a spectrum with a *.spec file as follows. 


beam num=1 x.origin=0 y.origin=$lightY angle=90 back.refl 
power.file=1300BB.spec wavel.start=0.207 wavel.end=5 


wavel.num=1000 
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The *.spec file gives the spectral radiance at various wavelengths in the region of 
interest. This command tells ATLAS to read in the data from the file ‘1300BB.spec’ and 
to solve a cell using 1000 evenly spaced rays in the interval from 0.207 um—5 um. The 


command also defines the origin of the beams and the incident angle. 


9. Solving 

The cell and its environment are now compktely defined in ATLAS and the pro- 
gram can now use this model to solve. The equations that are to be solved depend on 
what type of output is desired. This thesis concentrates on current-voltage relationships 
and photogeneration rates. By providing ATLAS with a number of current levels, it can 


solve for the associated voltage levels in order to produce an J-V curve as follows. 
solve ianode=0 b1=$I 
solve ianode=-$il b1=$I 
solve ianode=-$i2 b1=$I 


solve ianode=-$i3 b1=$I 








solve ianode=-$i4 b1=$I 


If photoge neration rates are desired, the user can tell ATLAS to solve for the pho- 


togeneration with discrete wavelengths of light incident on the cell as follows. 
solve b1=0.1 lambda=0.3 
solve b1=0.1 lambda=0.4 
solve b1=0.1 lambda=0.5 


solve b1=0.1 lambda=0.6 


The output of the data generated from solving these equation is stored in a *.log 


file which can be viewed graphically using TonyPlot. 
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10. Source Code 
The sample code given in this section was taken from a model for a GaSb cell. 


The full source code for both this cell and an InGaAs cell can be found in Appendix A. 


C. MATLAB 


The Silvaco software package is capable of performing all of the actions neces- 
sary to model a cell and graph the output, but it was often useful to use Matlab programs 


to supplement Silvaco for some operations. 


1. Iterations 

While ATLAS is capable of solving many very complex equations, it does lack 
some very basic programming capabilities such as conditional statements and loops that 
are necessary for performing iterations. Several programs (found in Appendix B) were 
written in Matlab to execute these functions. The program ‘atlasrun.m’ was used to run 
the iterations. In each iteration, ‘atlasrun.m’ would run an input deck with ATLAS, read 
and analyze the output data, iterate a parameter in the input deck, and rerun ATLAS. The 
program ‘ivmaxp.m’ solved for open-circuit voltage, short circuit current, maximum 
power, fill factor, and efficiency for each iteration and recorded this data in an output file. 


‘filerw.m’ was used to iterate a specified parameter in the input deck after each run. 


2; Plotting 

The TonyPlot program provided by Silvaco is capable of producing plots of all 
data that is generated by ATLAS, but it often more convenient to use Matlab to plot the 
data. In addition to running the iteration scheme, ‘atlasrun.m’ would also update an effi- 
ciency plot after each run of ATLAS. Since most optimization schemes took several 
hours to run, this feature provided up to date visual feedback on the progress of the rou- 
tine. The program ‘filerw.m’ was also used to plot /-V curves with a marked maximum 


power point on each curve. 


All of these programs can be found in full in Appendix B. 
a3 


This chapter covered the simulation software used in this thesis. Silvaco and Mat- 
lab were both utilized extensively for building and analyzing TPV cell models. The next 


chapter will cover the actual models built along with their resulting outputs. 
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VI. MODELS AND RESULTS 


This chapter discusses the models that were built in ATLAS. In each case, a sam- 
ple model was built and the output of this model was compared to accepted values. Once 
the model was verified, it was optimized for both AMO and the spectrum from a 1300K 


blackbody using an iteration routine written in Matlab. 


A. MODEL CONFIRMATION 


Before a TPV cell could be optimized for a blackbody spectrum, it was first nec- 
essary to confirm that the model being used was accurate. Although there was no data 
available on cells under a 1300K blackbody spectrum, there was available data on many 
different cells under the AMO spectrum. In addition, using Equation 2.2 the frequency 
response of a cell could be estimated. The assumption is that if the model can accurately 
predict the output of a cell under AMO and has the proper frequency response, it should 


also accurately predict the output under a different spectrum. 


1. Gallium Arsenide 

Gallium Arsenide (GaAs), while not suitable for TPV applications, is one of the 
most widely used materials for solar cells and therefore the maximum current-voltage 
relationship for this cell is very well documented. A 1 cm’ GaAs cell should provide ap- 
proximately 7, =31 mA and V,.=1 V [Ref. 3]. The cell in Figure 27 was built and 
tested in ATLAS. ATLAS was found to do a poor job of modeling antireflective coat- 
ings, so a 5% loss was assumed in all models. With present AR coating technology this 


is a very reasonable assumption. 
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Figure 27. GaAs Solar Cell 


Figure 28 shows the IV curve from this cell. 
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Figure 28. GaAs Cell Results 
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From Equation 2.2, the spectral response of a GaAs cell should approach zero at a 
wavelength of approximately 0.89 um. Figure 29 shows the calculated spectral response 
with a peak at 0.87 um compared to the normalized AMO spectrum. As expected, this 
graph does approach zero at 0.89 um. It is also important to note that the spectral re- 


sponse is closely correlated to the spectrum, thus the large power output. 
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Figure 29. Spectral Response of GaAs 


It is important that this graph go to zero in the correct place. This may not affect 
some cells under AMO because only a relatively small amount of energy is in the spec- 
trum above where the response goes to zero, but the spectrum from a 1300K blackbody, 
however, contains most of its energy at wavelengths longer than the cell should be able to 
respond to. If the spectral response does not go to zero, the model will allow the cell to 


respond to these higher wavelengths and give incorrect results. 


This GaAs cell model demonstrates that ATLAS is capable of providing a very 


accurate model of a solar cell. 


of 


2. Gallium Antimonide 

The next cell built was a gallium antimonide (GaSb) cell. This material is cur- 
rently one of the best candidates for replacing the thermocouples currently used in 
RTG’s. It was difficult to find data to compare with the results of this model. Most pub- 
lications in this area measure these cells under concentrators and skew the efficiency cal 
culations to appear more favorable. One documented cell claimed to achieve approxi- 


mately V, =33 mV and J,, =34.7 mA under AM1.5 froma 1 cm’ cell [Ref. 24]. This 


cell model seen in Figure 30, measured under AMO produced the output shown in Figure 
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Figure 30. GaSb Solar Cell 
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Figure 31. GaSb Solar Cell /-V Characteristics 


In addition to testing this cell under AMO, the incident spectrum was altered to 
model the spectrum from a concentrator at 68 and 70 suns concentration [Ref. 25]. The 


results from these runs very closely correlated with published data on this cell. 


From Equation 2.2 the spectral response of a GaSb cell should approach zero at a 


wavelength of 1.7 um. Figure 32 shows the spectral response of a GaSb cell compared 


to the normalized AMO spectrum. It can be seen from the figure that this response peaks 


at A =1.65 um, and, as expected, it approaches zero at the correct wavelength. 
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Spectral Response vs. AMO 
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Figure 32. Spectral Response of GaSb 


3 Indium Gallium Arsenide 

The final cell built was an indium gallium arsenide. The bandgap of this material 
is related to the molar concentration of gallium in the compound. By varying this con- 
centration, the bandgap can be varied from 0.35 eV to 1.42 eV [Ref. 26]. At an In/Ga 
concentration of 0.72/0.28 InGaAs has an energy bandgap of 0.55 eV [Ref. 14]. It is this 


cell that was modeled here. 

Unlike GaSb, there was a great deal of published data on InGaAs cells under 
AMO and other spectra as well. Unfortunately, because research on this material is in the 
early stages, there is still a great deal of inconsistency in published results. The recorded 
short circuit current data varied greatly from 30mA/cm? to 60mA/cm’ and recorded 
open-circuit voltage data varied from 0.2 V to 0.4 V, so it was difficult to determine the 


accuracy of this model [Ref. 26]. The cell in Figure 33 was built and found to give re- 


sults that fell within this range. The actual output of this model can be seen in Figure 34. 
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Figure 33. InGaAs Solar Cell 
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Figure 34. InGaAs Solar Cell /-V Characteristics 


0.25 


The spectral response of this InGaAs cell should approach zero at a wavelength of 


2.25 um according to Equation 2.2. Figure 35 shows the computed spectral response of 


this cell. The spectral response peaks at A =1.55 um and approaches zero at the correct 
wavelength. 
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Spectral Response vs. AMO 
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Figure 35. Spectral Response of InGaAs 


All of the cells tested provided an accurate model under the AMO spectrum and 


produced the expected spectral responses. 


B. THERMOPHOTOVOLTAIC CELLS 


The previous chapter established the validity of the cell models. The next step 
was to take these models and optimize them for the spectrum of a 1300K blackbody. In 
order to do this, the Matlab programs previously discussed were used to run iteration rou- 


tines to find the optimum thicknesses and doping levels for each layer of a cell. 


1. Gallium Antimonide 


Figure 36 below shows the optimum GaSb cell found. 
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Figure 36. GaSb TPV Cell 


Figure 37 shows the output from this cell. 
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Figure 37. GaSb TPV Cell /-V Characteristics 


The photogeneration rate was also measured from this cell. Figure 38 shows the 


photogeneration rate at incrementally increasing wavelengths from 4 =0.2 um—2.0 um. 
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Figure 38. GaSb Photogeneration Rate 


In this picture, brighter colors represent higher photogeneration rates. It can be 
seen that longer wavelengths of light are absorbed deeper in the cell and cause more pho- 
togeneration to occur. At very short and very long wavelengths of light, however, there 


is almost no photogeneration. 


2s Indium Gallium Arsenide 


Figure 39 shows the optimum InGaAs cell found. 
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Figure 39. InGaAs TPV Cell 


Figure 40 shows the output of this cell when the spectrum of 1300K blackbody 


was applied to it. 
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Figure 40. InGaAs TPV Cell /-V Characteristics 
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The results from this cell were compared to published data on a similar InGaAs 


cell under a 1225K blackbody spectrum. 


The photogeneration rate in this cell was measured in the same manner as the 
GaSb cell. Figure 41 shows the photogeneration for this cell with incrementally increas- 


ing wavelengths of light from 4 =0.2 um—2.3 um. 
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Figure 41. InGaAs Photogeneration Rate 


Again it can be seen that longer wavelengths of light are absorbed deeper in the 
cell and produce more photogeneration until the incident light no longer has enough en- 
ergy to induce photogeneration. From these photogeneration figures and the previous 
spectral response graphs, it can be deduced that cells generally absorb light most effi- 


ciently that has energy only slightly higher than the energy bandgap of the material. 


This chapter compared simulated cell results under the AMO spectrum to meas- 
ured data in order to verify the models. The cells were then optimized for the spectrum 


from a 1300K blackbody and the results of these cells were presented. The next chapter 
66 


discusses conclusions reached from this thesis and recommendations for future work in 


this area. 
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Vil. CONCLUSIONS AND RECOMMENDATIONS 


This thesis developed and demonstrated reliable models of GaSb and InGaAs 
thermophotovoltaic cells for both the AMO spectrum and the spectrum of 1300K black- 
body. In addition to the Silvaco models, several programs were written to enable Matlab 
to communicate with Silvaco for running simulations and analyzing data. With the de- 
velopment of models for other cell materials and possible multijunction configurations, 
Silvaco’s ATLAS software could prove to be an invaluable tool for TPV industry and for 


the cell designers. 


This thesis assumed the heat source to be a 1300K blackbody, but with very little 
effort any other spectrum could be applied to these cells. The Matlab program developed 
here can produce the spectrum from any temperature blackbody over any range of wave- 
lengths, calculating the radiance at each wavelength along with the total power under the 


curve. The optimization routines would run exactly the same way. 


The models built here could easily be configured to model other cell materials as 
well. The most difficult part of building an accurate cell model is providing accurate ma- 
terial properties. These properties are difficult to find for some of the newer, more exotic 
ternary and quaternary materials. Optical parameters, one of the most important material 
properties for building an accurate model, can often only be found for wavelengths of 
light up to 1 or 2 um but, for TPV applications, this data is necessary for much longer 
wavelengths of light. Silvaco, while a powerful simulation tool, cannot accurately solve 


all the equations to model a cell without the proper variables. 


Once initial cell models are built, however, they can provide a cell designer in- 
stant feedback on the impact of changes made in any given design. With the lack of ac- 
curate TPV cell models, the current design process consists of building and testing many 
different cells in order to find the most efficient design. The method presented in this 
thesis allows a cell designer the option of using any configuration of any definable mate- 
rials with any conceivable doping schemes. Other ongoing research is even adding the 
capability of predicting the impact of radiation on a cell. No other model has shown the 
versatility and ease of design for TPV cells as the Silvaco software used here. 
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The next step in this research is to model a multijunction TPV cell. Recent re- 
search by Green [Ref. 27] and Michalopoulos [Ref. 28] has proven that Atlas is capable 
of modeling and optimizing multijunction solar cells. A multijunction configuration 
would have the potential to greatly improve the heat to electrical conversion efficiency of 


the TPV cell as well. 
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APPENDIX A: ATLAS CODE 


This appendix contains the Silvaco input decks that were used to build the models 


in this thesis for GaAs, GaSb, and InGaAs. 


A. GALLIUM ARSENIDE 


1. I-V Curve 


go atlas 


PE HE HE EE EEE FE EE EEE EE HE EEE EEE HE EE EH 


Parameters 


Cell Width (um) 
set cellwidth=500 


Cap/Contact (%) 
set contactpercent=8 


Standard # of Divisions 
set divs=10 


Contact 
set contactthickness=.1 


Cap 
t capthickness=.3 
set capdoping=3e18 














Cathode 
set cathodethickness=.2 
set cathodedoping=1.5e18 
Anode 


anodethickness=200 
anodedoping=7el7 





ct ct 





Back Surface Field (BSF) 
t bsfthickness=.5 
set bsfdoping=8e18 














HEE EEE EE EE EE EE EE EE EEE EEE EEE 


t cellwdiv=Scellwidth/Sdivs 
set 3D=100e6 /$cellwidth 

t capwidth=.01*Scontactpercent*Scellwidth 
set capwdiv=Scapwidth/Sdivs 








set bsfL=0 
set bsfH=SbsfL-Sbsfthickness 
set bsfdivs=S$bsfthickness/S$divs 





set anoL=SbsfH 
set anoH=SanoL-Sanodethickness 
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set anodivs=Sanodethickness/S$divs 
set catL=SanoH 

set catH=ScatL-Scathodethickness 
set catdivs=$cathodethickness/S$divs 
set capL=$catH 

set capH=ScapL-Scapthickness 

set capdivs=Scapthickness/S$divs 

set contL=ScapH 

set contH=ScontL-Scontactthickness 
set contdivs=Scontactthickness/Sdivs 
set contcenter=$cellwidth/2 

set contRR=Scontcenter+tScapwidth/2 
set contLL=Scontcenter-Scapwidth/2 
set lightY=ScontH-5 


mesh width=$3D 


#X—Mesh 
x.mesh 
x.mesh 
x.mesh 
x.mesh 


#Y—Mesh 
.mesh 
.mesh 
-mesh 
-mesh 
.-mesh 
-mesh 


KKK KKK 


#Region 
##Conta 
region 
y.max=$ 
##Caps 
region 


region 
region 
y.max=$ 


#Cathode 


region 
#Anode 
region 
#BSF 

region 


##ELect 
lectro 





loc=0 spac=Scellwdiv 
loc=ScontLL spac=$capwdiv 
loc=S$contRR spac=$capwdiv 
loc=$cellwidth spac=$cellwdiv 


loc=ScontH spac=S$contdivs 
loc=ScapH spac=Scapdivs 
loc=ScatH spac=Scatdivs 
loc=SanoH spac=Sanodivs 
loc=SbsfH spac=Sbsfdivs 
loc=SbsfL spac=Sbsfdivs 


Ss 
chs 


num=1 material=Gold x.min=ScontLL x.max=ScontRR y.min=ScontH 


contL 


num=2 material=GaAs x.min=ScontLL x.max=ScontRR y.min=ScapH y.max=ScapL 
##Vacuums 


num=3 material=Vacuum x.min=0 


x.max=ScontLL y.min=S$contH y.max=ScapL 





num=4 material=vacuum x.min=$contRR x.max=Scellwidth y.min=$contH 


capL 


num=5 material=GaAs x.min=0 x. 


num=6 material=GaAs x.min=0 x. 





num=7 material=GaAs x.min=0 x. 


rodes 
de name=cathode material=Gold 





y.max=$ 
lectro 


contL 
de name=anode material=Gold x. 





y.max=SbsfL 


max=Scellwidth y.min=ScatH y.max=ScatL 
max=Scellwidth y.min=SanoH y.max=SanoL 


max=Scellwidth y.min=SbsfH y.max=SbsfL 


x.min=ScontLL x.max=ScontRR y.min=ScontH 


min=0 x.max=Scellwidth y.min=SbsfL 
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Doping 

Cap 
doping uniform region=2 n.type conc=Scapdoping 
Cathode 
doping uniform region=5 n.type conc=Scathodedoping 
Anode 
doping uniform region=6 p.type conc=Sanodedoping 
BSF 
doping uniform region=7 p.type conc=Sbsfdoping 








#Material properties 
material COPT=7e-10 
# Vacuum 
material material=Vacuum real.index=3.3 imag.index=0 
# GaAs 
material material=GaAs EG300=1.42 PERMITTIVITY=13.1 AFFINITY=4.07 
material material=GaAs NC300=4.7e17 NV300=7e18 

material material=GaAs index.file=GaAs.opt 

# Gold 
material material=Gold real.index=1.2 imag.index=1.8 
# Models 

models region=3 CONMOB 

models region=4 CONMOB 

models region=8 CONMOB 

models region=9 CONMOB 

models region=10 CONMOB 

models OPTR print 




















#Light beams 
beam num=1 x.origin=0 y.origin=SlightY angle=90 back.refl 
power.file=AMOsilv.spec wavel.start=0.21 wavel.end=4 wavel.num=50 





### SOLVING I-V CURVE 





set suns=1 
set I=Ssuns*.95 
set a=SI*.9 


solve init 
method gummel newton maxtraps=10 itlimit=25 


solve bl=Sa 


method newton maxtraps=10 itlimit=100 
solve b1=SI 


extract name="isc" max(curve(v."anode", i."cathode") ) 








set isc=Sisc*$3D 

set i1=Sisc/10 

set i2=$i1+$isc/10 
set 13=$i2+Sisc/10 
set i4=$13+$isc/10 
set i15=$i4+Sisc/10 
set i6=$i5+S$isc/20 
set i7=$16+$isc/20 
set i18=$i7+Sisc/20 
set i9=$i8+Sisc/20 
set i110=$i9+Sisc/20 
set i111=$i10+Sisc/40 
set i112=$i11+$isc/40 
set 113=$i112+Sisc/40 
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set 114=$i13+Sisc/40 
set 115=$i14+Sisc/40 
set i16=$115+Sisc/80 
set i17=$i116+$isc/80 
set i18=$117+Sisc/80 
set i19=$118+Sisc/80 
set i20=$119+Sisc/80 
set 121=$i20+Sisc/80 
set 122=$i21+Sisc/80 
set i23=$122+Sisc/80 
set 124=$i23+Sisc/80 
set i25=$124+Sisc/80 


contact name=anode current 
method newton maxtraps=4 itlimit=100 


log outfile=gaas.log 
struct outfile=gaas.str 


solve ianode=0 bi=SI 






































solve ianode=-$il b1=$I 
solve ianode=-$i2 b1=$I 
solve ianode=-$i3 b1=$I 
solve ianode=-$i4 b1=$I 
solve ianode=-$i5 b1=$I 
solve ianode=-$i6 b1=$I 
solve ianode=-$i7 b1=$I 
solve ianode=-$i8 b1=$I 
solve ianode=-$i9 b1=SI 
solve ianode=-$i10 b1=$I 
solve ianode=-$ill b1=$I 
solve ianode=-$i12 b1=$I 
solve ianode=-$i13 b1=$I 
solve ianode=-$i14 b1=$I 
solve ianode=-$i15 b1=$I 
solve ianode=-$il16 b1=$I 
solve ianode=-$il7 b1=$I 
solve ianode=-$i18 b1=$I 
solve ianode=-$i19 b1=SI 
solve ianode=-$i20 b1=$I 
solve ianode=-$i21 b1=$I 
solve ianode=-$i22 b1=$I 
solve ianode=-$i23 b1=$I 
solve ianode=-$i24 b1=$I 
solve ianode=-$i25 b1=$I 





log off 





tonyplot gaas.log -set onecon.set 
tonyplot gaas.str 


2. Frequency Response 


go atlas 


Hat aa at HE aE HE aE HE aE HE EE EE HE aE aE EE HE aE aE HE EE HE EE EE aE EE EE 
## Parameters 

Ft 

## Cell Width (um) 
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set cellwidth=500 











Cap/Contact (%) 


set contactpercent=8 


Standard # of Divisions 


set divs=10 























Contact 

set contactthickness=.1 
Cap 

set capthickness=.3 

set capdoping=3e18 
Cathode 

set cathodethickness=.2 

set cathodedoping=1.5e18 
Anode 

set anodethickness=200 

set anodedoping=7el7 


Back Surface Field (BSF) 
t bsfthickness=.5 
t bsfdoping=8el18 


HEE aE EE EE EE HE HE EEE EEE EEE EEE 



































set cellwdiv=Scellwidth/$divs 

set 3D=100e6 /$cellwidth 

set capwidth=.01*S$contactpercent*$cellwidth 
set capwdiv=Scapwidth/Sdivs 

set bsfL=0 

set bsfH=SbsfL-Sbsfthickness 

set bsfdivs=Sbsfthickness/Sdivs 

set anoL=SbsfH 

set anoH=SanoL-Sanodethickness 

set anodivs=Sanodethickness/S$divs 
set catL=SanoH 

set catH=ScatL-Scathodethickness 
set catdivs=$cathodethickness/S$divs 
set capL=$catH 

set capH=ScapL-Scapthickness 

set capdivs=$capthickness/S$divs 

set contL=ScapH 

set contH=ScontL-Scontactthickness 
set contdivs=Scontactthickness/Sdivs 
set contcenter=S$cellwidth/2 

set contRR=Scontcenter+Scapwidth/2 
set contLL=Scontcenter-Scapwidth/2 
set lightY=ScontH-5 





mesh width=$3D 
#X—Mesh 
x.mesh loc=0 spac=Scellwdiv 
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x.mesh loc=ScontLL spac=Scapwdiv 
x.mesh loc=ScontRR spac=Scapwdiv 
x.mesh loc=Scellwidth spac=$cellwdiv 


#Y—Mesh 


-mesh 
.mesh 
-mesh 
-mesh 
.-mesh 
.-mesh 


KKKK KK 


loc=ScontH spac=$contdivs 


loc=ScapH 
loc=S$catH 
loc=SanoH 
loc=SbsfH 
loc=SbsfL 


spac=Scapdivs 
spac=Scatdivs 
spac=Sanodivs 
spac=Sbsfdivs 
spac=Sbsfdivs 


#Regions 
##Contacts 
region num=1 
y.max=ScontL 
##Caps 
region num=2 
##Vacuums 
region num=3 
region num=4 
y.max=ScapL 





ma 








material=Gold x. 


material=GaAs x. 


material=Vacuum 
terial=vacuum 


-min=0O x. 


-min=0O x. 


-min=0O x. 


rial=Gold 





min=ScontLL x.max=$contRR y.min=S$contH 





min=ScontLL x.max=$contRR y.min=ScapH y.max=ScapL 





x.min=0 x.max=ScontLL y.min=ScontH y.max=ScapL 
x.min=ScontRR x.max=Scellwidth y.min=ScontH 





max=Scellwidth y.min=ScatH y.max=ScatL 
max=Scellwidth y.min=SanoH y.max=SanoL 


max=Scellwidth y.min=SbsfH y.max=SbsfL 





#Cathode 
region num=5 material=GaAs x 
#Anode 
region num=6 material=GaAs x 
#BSF 
region num=7 material=GaAs x 
##Electrodes 
lectrode name=cathode mat 
y.max=ScontL 
lectrode name=anode material=Gold x. 


*x.min=ScontLL x.max=ScontRR y.min=ScontH 





y-max=SbsfL 


Doping 

Cap 
doping uniform 
Cathode 
doping uniform 
Anode 
doping uniform 
BSF 
doping uniform 








region=2 n 
region=5 
region=6 p 


region=7 p 


#Material properties 


material 
# Vacuum 
material 
# GaAs 


COPT=7e-10 


.type 


-type 


.type 


.type 


min=0 x.max=Scellwidth y.min=SbsfL 


conc=Scapdoping 
conc=Scathodedoping 
conc=Sanodedoping 


conc=Sbsfdoping 


material=Vacuum real.index=3.3 imag.index=0 


ma 
ma 
ma 
# 

ma 


# 


terial 
terial 
terial 
Gold 
terial 
odels 








material=GaAs 
material=GaAs 
material=GaAs 


material=Gold 








EG300=1.42 PERMITTIVITY=13.1 AFFINITY=4.07 
C300=4.7e17 NV300=7e18 


index. file=GaAs.opt 





models region=3 CONMOB 


real.index=1.2 imag.index=1.8 
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models 
models 
models 
models 
models 


#Light 


region=4 CONMOB 
region=8 CONMOB 
region=9 CONMOB 
region=10 CONMOB 
OPTR print 


beams 


beam num=1 x.origin=0 y.origin=SlightY angle=90 back.refl 
power.file=AMOsilv.spec wavel.start=0.21 wavel.end=4 wavel.num=50 


### SO 
Seb. Su 
set I= 


set. a= 


solve 


LVING I-V CURVE 





ns=1 
Ssuns*.95 
SI*.9 


init 





method gummel newton maxtraps=10 itlimit=25 


solve 


b1l=Sa 


method newton maxtraps=10 itlimit=100 


solve 


extract name="isc" max(curve(v."anode", 


output 
solve 


log ou 


Solve 
solve 
solve 
solve 
Solve 
solve 
solve 
Solve 
Solve 
solve 
solve 
solve 
solve 





b1=SI 





con.band val.band 
init 
tfile=freq.log 
b1=0.1 lambda=0.22 
b1=0.1 lambda=0.225 
b1=0.1 lambda=0.23 
b1=0.1 lambda=0.235 
b1=0.1 lambda=0.24 
b1=0.1 lambda=0.245 
b1=0.1 lambda=0.25 
b1=0.1 lambda=0.255 
b1=0.1 lambda=0.26 
b1=0.1 lambda=0.265 
b1=0.1 lambda=0.27 
b1=0.1 lambda=0.275 
b1=0.1 lambda=0.28 
b1=0.1 lambda=0.285 
b1=0.1 lambda=0.29 
b1=0.1 lambda=0.295 
b1=0.1 lambda=0. 30 
b1=0.1 lambda=0.305 
b1=0.1 lambda=0.31 
b1=0.1 lambda=0.315 
b1=0.1 lambda=0.32 
b1=0.1 lLambda=0.325 
b1=0.1 lambda=0.33 
b1=0.1 lLambda=0.335 
b1=0.1 lambda=0.34 
b1=0.1 lambda=0.345 
b1=0.1 lambda=0.35 
b1=0.1 lLambda=0.355 
b1=0.1 lambda=0.36 
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i."cathode") ) 


sol 
sol 
sol 
sol 
sol 
sol 
Sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 





ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
Ve 
ve 
ve 
ve 
ve 
ve 
Ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 





lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 








lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
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sol 
sol 
sol 
sot 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 


so 


sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 


So 


sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 


So 


sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 
sol 


ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
lve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
lve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
lve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 
ve 





b1=0 
b1=0 
bl= 
bl1l= 
bl= 
b1=0 
b1=0 
b1=0 
bl= 
b1=0 
b1=0 
b1=0 
b1=0 
b1=0 
b1=0 
b1=0 
b1=0 
b1=0 
bl= 
b1=0 
b1=0 
b1=0 
bl= 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0 
bl 
b1=0. 
b1=0. 
b1=0. 
bl 
b1=0. 
bl= 
b1=0. 
bl= 
b1=0. 
bl= 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
b1=0. 
bl= 
b1=0. 
b1=0. 
b1=0 
b1=0 





PRP RPRPPPRPPRPRBPEPPEPHEPRPEPHPHEHP EEE 


PRP ERP RPP PPP PPE EY 


PPRPPPRPP Pt 


bh +t 

















lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0 
lambda=0. 
Lambda=0. 
lambda=0. 
ambda=0. 
lambda=0. 
ambda=0 
ambda=0 
ambda=0. 
lambda=0. 
ambda=0. 
lambda=0. 
ambda=0. 
lambda=0. 
ambda=0. 
lambda=0. 
ambda=0. 
ambda=0. 
ambda=0. 
lambda=0. 
ambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
ambda=0. 
lambda=0. 
ambda=0. 
ambda=0. 
ambda=0. 
lambda=0. 
ambda=0. 
lambda=0. 
ambda=0. 
lambda=0. 
ambda=0. 
lambda=0. 
ambda=0 
ambda=0 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
lambda=0. 
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solve b1=0.1 lambda=0.995 
solve b1=0.1 lambda=1.00 
log off 


B. GALLIUM ANTIMONIDE 


1. I-V Curve 


go atlas 


HHPEEEEEEEE EH HHH HH EE HE EE EE HE HE EE HE EE HE HE HE HE HE EE HE EE EE 
## Parameters 

## 
## Cell Width (um) 
set cellwidth=500 
## 
## Cap/Contact (%) 
set contactpercent=8 
## 
## Standard # of Divisions 
set divs=20 











## 

## Contact 

set contactthickness=.1 

# # 

## Cap (there is no cap in this model) 


set capthickness=.1 
set capdoping=3el18 
## 
## Cathode (N-type) 

set cathodethickness=10 

set cathodedoping=1.000000e+020 
## 
## Anode (P-type) (.2-.5um) 
set anodethickness=.37 

set anodedoping=9.8el19 

## 
## Back Surface Field (BSF) 

set bsfthickness=15.000000 

set bsfdoping=1.000000e+020 

## 

Ha a EH a a aE HE EE a aE EH EE a EH EE EE aE EE EE EE EE EE EE EH HE 














set cellwdiv=Scellwidth/S$divs 

set 3D=100e6 /Scellwidth 

set capwidth=.01*Scontactpercent*$cellwidth 
set capwdiv=Scapwidth/S$divs 

set bsfL=0 

set bsfH=SbsfL-Sbsfthickness 

set bsfdivs=Sbsfthickness/Sdivs 

set catL=SbsfH 

set catH=S$catL-S$cathodethickness 
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set catdivs=Scathodethickness/$divs 
set anoL=S$catH 

set anoH=SanoL-Sanodethickness 

set anodivs=Sanodethickness/Sdivs 
set capL=SanoH 

set capH=ScapL-Scapthickness 

set capdivs=Scapthickness/S$divs 

set contL=S$capH 

set contH=S$contL-S$contactthickness 
set contdivs=Scontactthickness/Sdivs 
set contcenter=$cellwidth/2 

set contRR=Scontcenter+Scapwidth/2 
set contLL=Scontcenter-Scapwidth/2 
#Sets light source 200um above surface (try different heights) 
set lightY=ScontH-200 














mesh width=$3D 

#X—-Mesh 

x.mesh loc=0 spac=Scellwdiv 
x.mesh loc=S$contLL spac=$capwdiv 
x.mesh loc=$contRR spac=S$capwdiv 
x.mesh loc=$cellwidth spac=$cellwdiv 








#Y—Mesh 

y.-mesh loc=$contH spac=Scontdivs 
y.mesh loc=$capH spac=Scapdivs 
y.mesh loc=$catH spac=Scatdivs 
y.mesh loc=$anoH spac=Sanodivs 
y.mesh loc=$bsfH spac=Sbsfdivs 
y.mesh loc=S$bsfL spac=Sbsfdivs 








#Regions 

##Contacts 
region num=1 material=Gold x.min=$contLL x.max=ScontRR y.min=ScontH 
y.max=ScontL 
##Cap 
region num=2 material=GaSb x.min=ScontLL x.max=ScontRR y.min=ScapH 
y.max=ScapL 
##Vacuums 
region num=3 material=Vacuum x.min=0 x.max=ScontLL y.min=ScontH 
y.max=ScapL 




















region num= aterial=vacuum x.min=ScontRR x.max=S$cellwidth 
y.min=ScontH y.max=ScapL 


#Anode 
region num=5 material=GaSb x.min=0 x.max=Scellwidth y.min=SanoH 
y.max=SanoL 
#Cathode 
region num=6 material=GaSb x.min=0 x.max=Scellwidth y.min=ScatH 
y.max=ScatL 
#BSF 
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region num=7 material=GaSb x.min=0 x.max=$cellwidth y.min=$bsfH 


y.max=SbsfL 


Hi! 


(T] 


lectrodes 





lectrode name=anod 


material=Gold x.min=$contLL x.max=S$contRR 








y.-min=ScontH y.max=ScontL 

















y.min=SbsfL y.max=SbsfL 






















































































#Doping (In this model 
eral) 

# Cap 

doping uniform region=2 
# Anode 

doping uniform region=5 
# Cathode 

doping uniform region=6 
# BSF 

doping uniform region=7 
#Material properties 

# Vacuu 

naterial aterial=Vacuu 
# GaSb 

material material=GaSb 
material material=GaSb 
naterial aterial=GaSb 
materia aterial=GaSb 
# Gold 

naterial aterial=Gold 
# Models 

models OPTR print 
#Light beams 








lectrode name=cathode material=Gold x.min=0 x.max=Scellwidth 





there are 


no specified doping levels, try sev- 





nc=Scapdoping 








p.type co 
p.type co 
n.type co 
n.type co 
COPT=5e-10 
EG300=.726 





index.file 


real.index 


nc=Sanodedoping 


nc=Scathodedoping 





nc=Sbsfdoping 


real.index=3.8 imag.index=0 


PERMITTIVITY=14.4 AFFINITY=4.06 





NC300=2.1le17 NV300=1.8e19 
=GaSb.opt 


=1.2 imag.index=1.8 


beam num=1 x.origin=0 y.origin=S$lightY angle=90 back.refl 
power.file=1300BB.spec wavel.start=0.207 wavel.end=5 wavel.num=1000 


### SOLVING I-V CURV 


set suns=1 
set I=Ssuns*.95 
set a=SI*.9 


solve init 





solve bl=Sa 





Ld 





1 
method gummel newton maxtraps=10 itlimit=25 
1 


method newton maxtraps=10 itlimit=100 


solve bl1=SI 








extract name="isc" 





set isc=S$isc*S$3D 


max (curve(v."anode", i."cathode") ) 
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set il=Sisc/10 

set i2=$11+$isc/10 
set i3=$i2+Sisc/10 
set i4=$i3+S$isc/10 
set i5=$14+Sisc/10 
set i16=$15+Sisc/20 
set i7=$16+Sisc/20 
set i8=$17+Sisc/20 
set 19=$i8+Sisc/20 
set 0=$19+Sisc/20 
set 111=$110+Sisc/40 
set 112=$i11+Sisc/40 
set i13=$112+Sisc/40 
set i114=$113+Sisc/40 
set 115=$114+Sisc/40 
set i16=$115+Sisc/80 
set i17=$116+Sisc/80 
set 118=$i17+Sisc/80 
set i19=$i18+Sisc/80 
set i20=$119+Sisc/80 
set 121=$i120+Sisc/80 
set 122=$121+Sisc/80 
set i123=$122+Sisc/80 
set 124=$i123+Sisc/80 
set i25=$124+Sisc/80 








contact name=anode current 





method newton maxtraps=4 itlimit=100 


log outfile=gasb.log 
struct outfile=gasb.str 





















































solve ianode=0 b1=SI 

solve ianode=-$il b1=SI 
solve ianode=-$i2 b1=SI 
solve ianode=-$i3 b1=SI 
solve ianode=-$i4 b1=SI 
solve ianode=-$i5 b1=SI 
solve ianode=-$i6 b1=SI 
solve ianode=-$i7 b1=SI 
solve ianode=-$i8 b1=SI 
solve ianode=-$i9 b1=SI 
solve ianode=-$i10 b1=SI 
solve ianode=-$ill b1=SI 
solve ianode=-$il2 b1=SI 
solve ianode=-$i13 b1=SI 
solve ianode=-$il14 b1=SI 
solve ianode=-$i15 b1=SI 
solve ianode=-$il16 b1=SI 
solve ianode=-$il7 b1=SI 
solve ianode=-$i18 b1=SI 
solve ianode=-$il9 b1=SI 
solve ianode=-$i20 b1=SI 
solve ianode=-$i21 b1=SI 
solve ianode=-$i22 b1=SI 
solve ianode=-$i23 b1=SI 
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solve ianode=-$i24 b1=SI 
solve ianode=-$i25 b1=SI 


log off 








tonyplot gasb.log -set onecon.set 
tonyplot gasb.str 





Ze 


go atlas 


Cell 


Cap 














Width 


Contact 
set contactthickness=.1 


Frequency Response 


HEE HE HE HE EHH HE EE EEE EEE EEE EEE EEE HE EE BE 


Parameters 


(um) 


set cellwidth=500 


Cap/Contact (%) 
set contactpercent=8 


Standard # of Divisions 
set divs=20 


(there is no cap in this model) 


set capthickness=.1 
set capdoping=3e18 


## Cathode (N-type) 


se 
se 
tt 


cv st 


cathodethickness=10 
cathodedoping=1.000000e+020 


## Anode (P-type) (.2-.5um) 


se 
se 
tt 


## Back Surf 
bsfthick 
bsfdopin 


se 
se 
tt 


Ge et 





HRT HE EH HE 


anodethi 
anodedoping=9.8el19 











ckness=.37 


ace Field (BSF) 
ness=15.000000 
g=1.000000e+020 


HEH ETE EEE HE THEE EEE EE EEE HEE EH 


cellwdiv=$cellwidth/$divs 








3D=50e6 /$cellwidth 


set capwidth=.01*Scontactpercent*$cellwidth 
set capwdiv=Scapwidth/Sdivs 























set bsfL=0 

set bsfH=SbsfL-Sbsfthickness 

set bsfdivs=S$bsfthickness/S$divs 

set catL=SbsfH 

set catH=ScatL-S$cathodethickness 
set catdivs=$cathodethickness/S$divs 
set anoL=S$catH 

set anoH=SanoL-Sanodethickness 

set anodivs=Sanodethickness/Sdivs 
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capL=SanoH 





n 
ct ct ct 


n 
oO 
GRockoGh (Ch ocr Ck 


cont L=S$capH 
contH=S$contL-Scontactthickness 

contdivs=Scontactthickness/Sdivs 
contcenter=Scellwidth/2 
contRR=Scontcenter+Scapwidth/2 
contLL=Scontcenter-Scapwidth/2 


capH=ScapL-Scapthickness 
capdivs=Scapthickness/S$divs 


#Sets light source 200um above surface 
set lightY=ScontH-200 





mesh width=$3D 


#X—Mesh 





x.mesh loc=0 spac=Scellwdiv 
x.mesh loc=ScontLL spac=Scapwdiv 
x.mesh loc=$contRR spac=Scapwdiv 
x.mesh loc=Scellwidth spac=$cellwdiv 


#Y—Mesh 
y.mesh 
-mesh 
-mesh 
-mesh 
-mesh 
-mesh 


KKK KK 


loc=ScontH spac=S$contdivs 


loc=ScapH 
loc=ScatH 
loc=SanoH 
loc=SbsfH 
loc=SbsfL 


#Regions 
##Contacts 


region num=1 


y.max=ScontL 


##Cap 


region num=2 ma 
##Vacuums 


region num=3 


region num=4 ma 
y.max=ScapL 


Anode 


BSF 











region num=5 ma 
Cathode 
region num=6 ma 


region num=7 ma 


#ELlectrodes 


terial 


terial 


terial 


terial 


material=Gold x. 


l=GaSb x. 


material=Vacuum 


L=vacuum 


l=GaSb x. 


l=GaSb x. 








terial 


l=GaSb x. 


spac=Scapdivs 
spac=Scatdivs 
spac=Sanodivs 
spac=Sbsfdivs 
spac=Sbsfdivs 


min=ScontL 


min=ScontL 


(try different heights) 


L x.max=ScontRR y.min=S$contH 


L x.max=ScontRR y.min=ScapH y.max=ScapL 


x.min=0 x.max=ScontLL y.min=ScontH y.max=ScapL 


x.min=$con 





tRR x.max=Scellwidth y.min=ScontH 


min=0 x.max=Scellwidth y.min=SanoH y.max=SanoL 


min=0 x.max=Scellwidth y.min=ScatH y.max=ScatL 


min=0 x.max=Scellwidth y.min=SbsfH y.max=SbsfL 





lectrod 


y.max=ScontL 


name=anod 


material=Gold x.mi 





lectrod 


y-max=SbsfL 


#Doping 


# Cap 


(In this model ther 





n=ScontLL x.max=ScontRR y.min=ScontH 


name=cathode material=Gold x.min=0 x.max=Scellwidth y.min=SbsfL 





are no sp 


cified doping levels, try several) 


doping uniform region=2 p.type conc=Scapdoping 


# Anode 


doping uniform region=5 p.type conc=Sanodedoping 
# Cathode 
doping uniform region=6 n.type conc=S$cathodedoping 
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# BSF 
doping uniform region=7 n.type conc=Sbsfdoping 


material 
# Vacuum 
material 
# GaSb 
material 
material 
material 
# Gold 
material 





# 





odels 


COPT=7e-10 
material=Vacuum real.index=3.3 imag.index=0 
material=GaSb EG300=.72 PERMITTIVITY=15.7 AFFINITY=4.06 


material=GaSb NC300=2.07846e17 NV300=1.81865e17 
material=GaSb index.file=GaAs.opt 











material=Gold real.index=1.2 imag.index=1.8 


models OPTR print 


#Light beams 
beam num=1 x.origin=0 y.origin=SlightY angle=90 back.refl 
power.file=1300BB.spec wavel.start=0.1 wavel.end=7 wavel.num=100 


set I=.1 





solve init 
solve b1=S$I 


extract name="isc" max(curve(v."anode", i."cathode") ) 


output con.band val.band 
solve init 
log outfile=freq.log 














solve b1=0.1 lambda=0.22 
solve b1=0.1 lambda=0.3 
solve b1=0.1 lambda=0.4 
solve b1=0.1 lambda=0.5 
solve b1=0.1 lambda=0.6 
solve b1=0.1 lambda=0.7 
solve b1=0.1 lambda=0.8 
solve b1=0 lambda=. 82 
solve bl= lambda=. 84 
solve b1=0 lambda=. 86 
solve bl= lambda=. 88 
solve b1=0 lambda=0.9 
solve b1=0 lambda=1.0 
solve b1=0 lambda=1.1 
solve b1=0 lambda=1.2 
solve bl= lambda=1.3 
solve b1=0 lambda=1.4 
solve b1=0 lambda=1.5 
solve b1=0 lambda=1.6 
solve b1=0 lambda=1.7 
solve b1=0 lambda=1.8 
solve bl= lambda=1.9 
solve bl= lambda=2 .0 
solve bl= lambda=2.1 
solve b1=0 lambda=2.2 
solve bl= lambda=2.3 
solve bl= lambda=2.4 
solve bl= l lambda=2.5 
solve bl= lambda=2.6 
solve b1=0.1 lambda=2.7 
solve bl= lambda=2 .8 
solve bl= l lambda=2.9 
solve bl= lambda=3.0 
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solve b1=0.1 lambda=3.1 
solve b1=0.1 lambda=3.2 
solve b1=0.1 lambda=3.3 
solve b1=0.1 lambda=3.4 
solve b1=0.1 lambda=3.5 
solve b1=0.1 lambda=3.6 
solve b1=0.1 lambda=3.7 
solve b1=0.1 lambda=3.8 
solve b1=0.1 lambda=3.9 
solve b1=0.1 lambda=4.0 
solve b1=0.1 lambda=4.1 
solve b1=0.1 lambda=4.2 
solve b1=0.1 lambda=4.3 
solve b1=0.1 lambda=4.4 
solve b1=0.1 lambda=4.5 
solve b1=0.1 lambda=4.6 
solve b1=0.1 lambda=4.7 
solve b1=0.1 lambda=4.8 
solve b1=0.1 lambda=4.9 
solve b1=0.1 lambda=5.0 
solve b1=0.1 lambda=5.1 
solve b1=0.1 lambda=5.2 
solve b1=0.1 lambda=5.3 
solve b1=0.1 lambda=5.4 
solve b1=0.1 lambda=5.5 
solve b1=0.1 lambda=5.6 
solve b1=0.1 lambda=5.7 
solve b1=0.1 lambda=5.8 
solve b1=0.1 lambda=5.9 
solve b1=0.1 lambda=5.0 
solve b1=0.1 lambda=6.1 
solve b1=0.1 lambda=6.2 
solve b1=0.1 lambda=6.3 
solve b1=0.1 lambda=6.4 
solve b1=0.1 lambda=6.5 
solve b1=0.1 lambda=6.6 
solve b1=0.1 lambda=6.7 
solve b1=0.1 lambda=6.8 
solve b1=0.1 lambda=6.9 
solve b1=0.1 lambda=7.0 
log off 


C. INDIUM GALLIUM ARSENIDE 


1. I-V Curve 


go atlas 


HHPEEEEEEE EEE HEE HE EE EE HE EE HE HE HE HEH HE EE HE HE HE HE HE HE EE EE EE 
## Parameters 

## 

## Cell Width (um) 

set cellwidth=500 

## 

## Cap/Contact (%) 
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set contactpercent=8 


## Standard # of Divisions 
set divs=20 


## Contact 
set contactthickness=.1 








# # 

## 

##Molar Concentration of In(1-x)Ga(x)As (x) 
set x=.28 

set perm=15.1-2.87*Sx+.67*Sx%*2 

set affin=4.9-.83*Sx 

set n=3.51-.21*Sx 

# # 


## Cap (there is no cap in this model) 


set capthickness=.1 
set capdoping=3el18 
## 


## Cathode (N-type) 





set cathodethickness=17 

set cathodedoping=1e20 

Ht 

## Anode (P-type) (.2-.5um) 
set anodethickness=2.22 

set anodedoping=1e20 

## 


## Back Surface Field (BSF) 





set bsfthickness=.1 
set bsfdoping=1e20 
# # 


He aE aE aE aE ae aE aE AE aE a aE aE aE AE a HE HE ea aE aE HE aE aE aE aE aaa aE aaa aE aE 























set cellwdiv=Scellwidth/S$divs 

set 3D=100e6 /Scellwidth 

set capwidth=.01*Scontactpercent*$cellwidth 
set capwdiv=Scapwidth/S$divs 

set bsfL= 

set bsfH=SbsfL-Sbsfthickness 

set bsfdivs=Sbsfthickness/Sdivs 

set catL=SbsfH 

set catH=ScatL-Scathodethickness 

set catdivs=Scathodethickness/S$divs 
set anoL=ScatH 

set anoH=SanoL-Sanodethickness 

set anodivs=Sanodethickness/Sdivs 
set capL=SanoH 

set capH=ScapL-Scapthickness 

set capdivs=Scapthickness/S$divs 

set contL=S$capH 

set contH=S$contL-S$contactthickness 
set contdivs=Scontactthickness/Sdivs 
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set contcenter=Scellwidth/2 
set contRR=Scontcenter+Scapwidth/2 
set contLL=Scontcenter-Scapwidth/2 


#Sets light source 200um above surface (try different heights) 
set lightY=ScontH-200 





mesh width=$3D 

#X-Mesh 

x.mesh loc=0 spac=Scellwdiv 
x.mesh loc=S$contLL spac=S$capwdiv 
x.mesh loc=S$contRR spac=S$capwdiv 
x.mesh loc=S$cellwidth spac=Scellwdiv 








#Y—-Mesh 

y.mesh loc=$contH spac=Scontdivs 
y.mesh loc=$capH spac=Scapdivs 
y.-mesh loc=$catH spac=Scatdivs 
y.mesh loc=S$anoH spac=Sanodivs 
y.mesh loc=$bsfH spac=Sbsfdivs 
y.mesh loc=$bsfL spac=Sbsfdivs 








#Regions 

##Contacts 
region num=1 material=Gold x.min=$contLL x.max=ScontRR y.min=ScontH 
y.max=ScontL 
##Cap 
region num=2 material=InGaAs x.min=ScontLL x.max=ScontRR y.min=ScapH 
y.max=ScapL 
##Vacuums 

region num=3 material=Vacuum x.min=0 x.max=ScontLL y.min=ScontH 
y.max=ScapL 
region num= aterial=vacuum x.min=$contRR x.max=Scellwidth 
y.min=ScontH y.max=ScapL 
#Anode 
region num=5 material=InGaAs x.min=0 x.max=$cellwidth y.min=SanoH 
y.max=SanoL 
#Cathode 
region num=6 material=InGaAs x.min=0 x.max=$cellwidth y.min=$catH 
y.max=ScatL 
#BSF 
region num=7 material=InGaAs x.min=0 x.max=$cellwidth y.min=SbsfH 
y.max=SbsfL 

















































































































##ELectrodes 

lectrode name=anode material=Gold x.min=ScontLL x.max=ScontRR 
y.-min=S$contH y.max=ScontL 

lectrode name=cathode material=Gold x.min=0 x.max=$cellwidth 
y.-min=SbsfL y.max=SbsfL 





























#Doping (In this model there are no specified doping levels, try sev- 
eral) 
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doping uniform region=2 p.type conc=Scapdoping 

# Anode 

doping uniform region=5 p.type conc=Sanodedoping 

# Cathode 

doping uniform region=6 n.type conc=S$cathodedoping 
# BSF 

doping uniform region=7 n.type conc=$bsfdoping 































































































#Material properties 

# Vacuum 

material material=Vacuum real.index=3.4 imag.index=0 
# InGaSb 

material material=InGaAs COPT=7e-10 

naterial material=InGaAs EG300=.55 PERMITTIVITY=Sperm AFFINITY=Saffin 
materia aterial=InGaAs NC300=1.15e17 NV300=8.12e18 
material material=InGaAs index.file=InGaAs.opt 

# Gold 

material aterial=Gold real.index=1.2 imag.index=1.8 
# Models 

models OPTR print 











#Light beams 
beam num=1 x.origin=0 y.origin=S$lightY angle=90 back.refl 
power. file=1300bb.spec wavel.start=0.21 wavel.end=7 wavel.num=50 














### SOLVING I-V CURVI 





Pl 


set suns=1 
set I=.95*Ssuns 
set a=SI*.75 


solve init 
method gummel newton maxtraps=10 itlimit=25 
solve bl=Sa 


method newton maxtraps=10 itlimit=100 
solve b1=SI 








extract name="isc" max(curve(v."anode", i."cathode") ) 












































set isc=$isc*$3D 

set il=Sisc/10 

set i2=$11+$isc/10 
set i3=$12+Sisc/10 
set i4=$13+$isc/10 
set i5=$14+Sisc/10 
set i16=$15+Sisc/20 
set i7=$16+$isc/20 
set i8=$17+Sisc/20 
set i19=$18+Sisc/20 
set 110=$i19+Sisc/20 
set 111=$110+Sisc/40 
set 112=$i111+Sisc/40 
set i13=$112+Sisc/40 


90 























set 114=$113+Sisc/40 
set 115=$i114+Sisc/40 
set i16=$115+Sisc/80 
set 117=$i16+Sisc/80 
set 118=$i17+Sisc/80 
set i19=$118+Sisc/80 
set i20=$119+Sisc/80 
set 121=$i120+Sisc/80 
set 122=$i121+Sisc/80 
set 123=$i122+Sisc/80 
set i24=$123+Sisc/80 
set 125=$i124+Sisc/80 


contact name=anode current 





method newton maxtraps=10 itlimit=100 


log outfile=ingaas.log 






























































struct outfile=ingaas.str 
solve ianode=0 b1=SI 
solve ianode=-$il b1=SI 
solve ianode=-$i2 b1=SI 
solve ianode=-$i3 b1=SI 
solve ianode=-$i4 b1=SI 
solve ianode=-$i5 b1=SI 
solve ianode=-$i6 b1=SI 
solve ianode=-$i7 b1=SI 
solve ianode=-$i8 b1=SI 
solve ianode=-$i9 b1=SI 
solve ianode=-$i10 b1=SI 
solve ianode=-$ill b1=SI 
solve ianode=-$il2 b1=SI 
solve ianode=-$i13 b1=SI 
solve ianode=-$i14 b1=SI 
solve ianode=-$i1l5 b1=SI 
solve ianode=-$il16 b1=SI 
solve ianode=-$il7 b1=SI 
solve ianode=-$i18 b1=SI 
solve ianode=-$il19 b1=SI 
solve ianode=-$i20 b1=SI 
solve ianode=-$i21 b1=SI 
solve ianode=-$i22 b1=SI 
solve ianode=-$i23 b1=SI 
solve ianode=-$i24 b1=SI 
solve ianode=-$i25 b1=SI 
log off 








tonyplot ingaas.log -set onecon.set 
tonyplot ingaas.str 
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go 


He aE aE aE aE ae aE aE aE AE a aE aE aE AE a HE HE aE aE aE HE a aE aE aE HE aE aa aE aE aE 


Ze Frequency Response 


atlas 

















## Parameters 

## 

## Cell Width (um) 

set cellwidth=500 

# # 

## Cap/Contact (%) 

set contactpercent=8 

## 

## Standard # of Divisions 

set divs=20 

## 

## Contact 

set contactthickness=.1 

## 

# # 

##Molar Concentration of In(1-x)Ga(x)As(x) 
set x=.28 

set perm=15.1-2.87*Sx+.67*Sx%*2 
set affin=4.9-.83*Sx 

set n=3.51-.21*Sx 

## Cap (there is no cap in this model) 
set capthickness=.1 

set capdoping=3el18 

# # 

## Cathode (N-type) 

set cathodethickness=17.000000 
set cathodedoping=1e20 

## 

## Anode (P-type) (.2-.5um) 
set anodethickness=2.22 

set anodedoping=1+20 

## 

## Back Surface Field (BSF) 
set bsfthickness=.1 

set bsfdoping=1e20 

## 


He aE aE aE aE aE aE aE aE aE a aE aE aE AE a HE HE aE aE aE HE aE aE aE aE HE aa a aE aE aE EEE 

















cellwdiv=Scellwidth/$divs 








3D=100e6 /Scellwidth 





, 9 ee aT 


capwdiv=Scapwidth/Sdivs 


bs 
bs 


L=0 
H=SbsfL-Sbsfthickness 
divs=S$bsfthickness/S$divs 


Fh Fh fh 


Cr tT er 


L=SbsfH 
H=ScatL-S$cathodethickness 
catdivs=S$cathodethickness/Sdivs 














GP “ce “ct 
Q 
re) 
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capwidth=.01*Scontactpercent*Scellwidth 























set anoL=S$catH 

set anoH=SanoL-Sanodethickness 

set anodivs=Sanodethickness/S$divs 
set capL=SanoH 

set capH=ScapL-Scapthickness 

set capdivs=Scapthickness/S$divs 

set contL=ScapH 

set contH=ScontL-Scontactthickness 
set contdivs=Scontactthickness/S$divs 
set contcenter=Scellwidth/2 

set contRR=Scontcenter+Scapwidth/2 
set contLL=Scontcenter-Scapwidth/2 
#Sets light source 200um above surface (try different heights) 
set lightY=ScontH-200 











mesh width=$3D 

#X—-Mesh 

x.mesh loc=0 spac=$cellwdiv 
x.mesh loc=S$contLL spac=S$capwdiv 
x.mesh loc=$contRR spac=S$capwdiv 
x.mesh loc=$cellwidth spac=Scellwdiv 








#Y—Mesh 

y.mesh loc=$contH spac=Scontdivs 
y.-mesh loc=$capH spac=Scapdivs 
y.mesh loc=$catH spac=Scatdivs 
y.mesh loc=$anoH spac=Sanodivs 
y.mesh loc=$bsfH spac=Sbsfdivs 
y.mesh loc=$bsfL spac=Sbsfdivs 








#Regions 

##Contacts 
region num=1l material=Gold x.m 
y.max=ScontL 
##Cap 
region num=2 material=InGaAs x.min=ScontLL x.max=ScontRR y.min=ScapH 
y.max=ScapL 
##Vacuums 

region num=3 material=Vacuum x.min=0 x.max=ScontLL y.min=ScontH 
y.max=ScapL 


Ju. 


n=ScontLL x.max=$contRR y.min=ScontH 























region num= aterial=vacuum x.min=$contRR x.max=Scellwidth 
y.min=ScontH y.max=ScapL 

#Anode 

region num=5 material=InGaAs x.min=0 x.max=$cellwidth y.min=SanoH 
y.max=SanoL 


region num=6 material=InGaAs x.min=0 x.max=$cellwidth y.min=$catH 



















































































region num=7 material=InGaAs x.min=0 x.max=$cellwidth y.min=$bsfH 















































































































































##Electrodes 
lectrode name=anode material=Gold x.min=ScontLL x.max=ScontRR 
y.-min=ScontH y.max=ScontL 
lectrode name=cathode material=Gold x.min=0 x.max=Scellwidth 
y.min=SbsfL y.max=SbsfL 
#Doping (In this model there are no specified doping levels, try sev- 
eral) 
# Cap 
doping uniform region=2 p.type conc=Scapdoping 
# Anode 
doping uniform region=5 p.type conc=Sanodedoping 
# Cathode 
doping uniform region=6 n.type conc=$cathodedoping 
# BSF 
doping uniform region=7 n.type conc=Sbsfdoping 
#Material properties 
# Vacuu 
material aterial=Vacuum real.index=3.4 imag.index=0 
# InGaSb 
materia aterial=InGaAs COPT=7e-10 
naterial material=InGaAs EG300=.55 PERMITTIVITY=Sperm AFFINITY=Saffin 
materia aterial=InGaAs NC300=1.15e17 NV300=8.12e18 
material material=InGaAs index.file=InGaAs.opt 
# Gold 
material material=Gold real.index=1.2 imag.index=1.8 
# Models 
models OPTR print 
#Light beams 
beam num=1 x.origin=0 y.origin=$lightY angle=90 back.refl 





power.file=1300bb.spec wavel.start=0.21 wavel.end=5 wavel.num=100 





### SOLVING I-V CURVE 





set I=.1 


solve init 
solve b1=SI 


extract name="isc" max(curve(v."anode", i."cathode") ) 








output con.band val.band 
solve init 

log outfile=freq.log 
solve b1=0.1 lambda=0.22 
solve bl= lambda=0.3 
solve bl= lambda=0.4 
solve bl= ambda=0.5 
solve bl= lambda=0.6 
solve bl= lambda=0.7 
solve bl= lambda=0.8 
solve bl= lambda=. 82 
solve bl= lambda=. 84 
solve bl= lambda=.86 
solve bl= lambda=. 88 
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APPENDIX B: MATLAB CODE 


This section contains the Matlab functions that were written and used for this the- 


SiS. 


A. ATLASARUN.M 


Satlasrun.m 


oe 


6This program runs iterations of Atlas. It uses ivmaxp.m to calculate 
open %Scircuit voltage, short circuit current, max power, fill factor, 
and efficiency. 

6The function filerw.m is called to iterate the ATLAS input deck. The 
function %Stime.m is called to calculate the time of each iteration in 
hours, minutes, and %seconds. 








%(c) 2004 by B.P. Davenport 








sspecifies the number of iterations to perform 
iterations=250; 

sspecifies the range of values to iterate between 
old=1.8;0l=old; 

final=5; 








step=(final-old) /iterations; 


k=0; 
totalruntime=0; 


for (a=ol:step:final) 
tic 
k=k+1; 
fidres=fopen('atres.txt','a') 
filerw('ingaas.in',old,a) 
delete done.log 
delete atlaslog.log 





sHere atlas is called by Matlab to run the specified input deck 
'D:\Silvaco\etc\GuiAppStarter.ex lib-dir-name deckbuild -exe-name 
Deckbld -run ingaas.in -outfile atlaslog.log 











sprintf ('Executing run %u\nStandby for results',k) 
xy=-l; 





while (xy==-1) 
xy=fopen('done.log'); 
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paus 
end 


e(3) 


SHere the *.log file created by ATLAS is read and open-circuit 


voltage, 
Sshort 
returned. 
[Voc,Isc 
data(k, 


:)=[a Voc Isc Pmax FF Eff]; 


circuit current, max power, fill factor and efficiency are 





,Pmax,FF,Eff]=ivmaxp('ingaas.log'); 





sThe previously calculated data is written to file. 











fprintf(fidres,'Sf Sf Sf Sf Sf Sf\n',data(k,:)); 
SAfter each iteration, all calculated data is printed to the screen 


° 


ciency is 





° 


Supdated 
ele 
data 











Sincluding the time it took to run the iteration. A plot of effi- 


after each iteration. 


thisruntime=toc;Tl=time (toc); 


totalrun 


time=totalruntimet+tthisruntime; T2=time (totalruntime) ; 


averageruntime=totalruntime/k;T3=time (averageruntime) ; 


estimate 
time=average 
sprintf ( 





run- 
runtime* (iterations+l1) ;T4=time (estimateruntime) ; 
"This run took %s\nTotal run time so far is %s\nAverage run 








time is %s\nEstimated total run time is %s',T1,T2,T3,T4) 


ld=a; 


ause (3) 





tO MO 





close('all') 


lot (data(:,1),data(:,6) ) 
itle('Efficiency') 





end 


igure (1) 
title('Voc"') 
igure (2) 
itle('Isc') 
Figure (3) 
title ('Pmax') 
Figure (4) 
title ('FF') 


Figure (5) 











title ('Eff') 





plot (data(:,1),data(:,2)) 
Spa ear aaa a 
plot (data(:,1),data(:,4)) 
plot (data(:,1),data(:,5)) 


plot (data(:,1),data(:,6)) 
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B. FILERW.M 


function filerw(file,old, new) 

ilerw(file, old, new) 

his program opens an infile, "file", and writes over the 

S"old" anode thickness with the "new" anodethickness 

File must be in ''! 

old and new are the old and new values of anodethickness 

This same function was also used to iterate other variables by replac- 


fo) 


ing %*‘anodethickness’ with the name of the variable to be iterated. 





4 
qd 














6(c) 2004 B.P. Davenport 


fidr=fopen(file,'r'); 

a=fscanf(fidr,'%c'); 

fFelose(fidr); 

fidw=fopen(file,'w'); 

a=strrep(a,sprintf ('anodethickness=%f',old),sprintf('anodethickness=%f' 
,new)); 

fwrite (fidw,a); 

fFclose (fidw); 





C. IVMAXP.M 





function [Voc,Isc,Pmax,FF,Eff]=ivmaxp (varargin) 
6[Voc,Isc,Pmax,FF,Eff]=ivmaxp (file, type) 

SFinds Voc, Isc, Pmax, FF, and Efficiency from an Atlas *.log file 
Sfile and type must be in '' 

sType is the doping of the top layer, n or p 








6(c) 2004 B.P. Davenport 
fF (nargin==1) 
file=char(varargin(1)); 
type='P'; 
elseif (nargin==2) 
file=char(varargin(1)) 
type=char (varargin(2) ) 











else 
"Type nus be either ENS Mansy Pe or oa 
Voc=NaN; 
Isc=NaN; 
Pmax=NaN; 
FF=NaN; 
return 

end 








if (type=='N'|type=='n') 

[I V]J=textread(file,'%S*s %*f S*f S*f S*E S*E S*E SE SE S*E 
S*f', 'headerlines',20); 
elseif (type=='P'|type=='p') 
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[V 
Sf', 'headerlines',20); 
else 

"Type mus be either 

Voc=NaN; 

Isc=NaN; 

Pmax=NaN; 

FF=NaN; 

return 
end 
P=1.*V; 
[a b]=max(P); 
Pmax=V (b) *I (b); 
Voc=max (V) ; 
Isc=max (1); 
FF=Pmax/ (Voc*Isc) ; 
Eff£=100*Pmax/ (.1353); 


UN 





x=1.1*Voc; 














y=1.1*Isc; 

plot (V,1I,V(b),I(b),'rh') 
xlabel('Voltage (V)') 
ylabel('Current (A)') 
title('I-V Curve') 
axis([0 x 0 yl) 

D. BLACKBODY.M 


function 


[Lx, Sx] =blackbody (spec) 

This program calculates 
sents AMO and 
G. 


I]=textread(file, '%S*s 


the 


S*f S*f Skt Skt 


uae "Pte or 


[Lx, Sx] =blackbody (spec) 


spectrum a 5800K 





he spectrum from any other temperatur 


body 








wavelength ‘sand radiance ar 
spectrum of the 
fo 


2004 B.P. Davenport 


spec=char (spec); 
Tsun=5800; 
T=1300; 





Sconstants definition 
h=6.6260693e-34; 
c=299792458; 
k=1.386505e-23; 


blackbody which repre- 


specified by ‘“T’. The 





returned for th 


desired spectrum. The 


Sblackbody source is written to a *.spec file formatted 
r ATLAS and the Sresults are plotted. 


scaling factor to account for the distance from the sun 





dilution=1/46200; 
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St 


he number of 





points=3000; 
L=linspace(le-18,11le-6,points) *10%6; 


Ssun=(2.*pi.*c.%2.*h) ./((L*10%-6) .*5.* (exp (h.*c./ ( (L*10%-6) 


1) 


Fi 


en 





)*le-10*dilution; 


)*le-10; 


xsun, ysun]=max(Ssun) ; 
(2.8 pin KE 2 4. * Aun J (at 10°36) "5: 2 Cexp (hee. /((b*L0%-6) 


xother, yother]=max(S); 


d=fopen('1300BB.spec','w') 





d 


Fclose (fid) 

















forintf (fid, 'Su\n',length(S) ) 
for (x=l:points) 
forintf(fid,'sf Sf\n',L(x),S(x)); 


(W/com*2/um) ') 


(W/cm*2/um) ') 


plot (L,Ssun,'b',L,S,'g"') 
title('Optical File') 
xlabel('Wavelength (um) ') 
ylabel('Power Density 
figure (2) 

plot (L,Ssun,'b') 
title('Sun Optical File') 
xlabel('Wavelength (um) ') 
ylabel('Power Density 
psun=0; 

pother=0; 





for (x=1: (points-1l) ) 


po 


psun=psun+ (L(x+1)—-L(x)) *Ssun (x) 
pother=pother+ (L(x+1)-L(x))*S (x); 


un 


cher 


cher 


fa 


venly spaced points to calculate 


otk. eT) ) = 


ps=sprintf('Total power is:\n%f Watts/cm*2\nThe peak is:\n%f 


Wa 


tts/com*2\nAt 


ps2=sprintf('Total 


Wa 








abel 


tts/com*2\nAt 





igure (3) 
plot (L,S,'g') 
itle('1300 Blackbody Radiance') 
("Wavelength (um) ') 





abe 


("Power Density 





Ce 











xE( 


Oy: ts pis) 


Figure (2) 
text (2,.1,ps2) 


Lx=L'; 
if (strcmp (spec, 'BB') |strcmp(spec, 'bb') ) 


Sx=S'; 


power is:\n%sf 


(W/com*2/um) ') 
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L=%f um',pother, xother,L(yother)); 


Watts/cm*2\nThe 


L=sf um',psun,xsun,L(ysun)); 


peak 


.*k.*Tsun) )—- 


is:\n%sf 


end 
if(strcemp (spec, 'Sun') |strcmp (spec, 'sun')) 
Sx=Ssun'; 


end 


E. FREQRESP.M 


function [L,photo]=freqresp (varargin) 


ol? 


6[L,photo]=freqresp (varargin) 
sThis program reads in data a specified file and returns and plots the 
frequency sresponse from a cell. 


° 


S(c) 2004 B.P. Davenport 
































SIf there is only one input, the cell is assumed to be P on N. 
if (nargin==1) 
file=char(varargin(1l)); 
type='P'; 
elseif (nargin==2) 
file=char(varargin(1l)); 
type=char(varargin(2)); 
elseif (nargin==3) 
file=char(varargin(1l)); 
type==char (varargin(2)); 
else 
"Type mus be either TONEN Naty "PA or tp! ' 
L=NaN; 
photo=NaN; 
return 
end 


This function can accept data in two different forms from ATLAS or 
from a file %containing wavelength radiance pairs. 

if (type=='N'|type=='n') 

L photo]=textread(file,'%*s %S*f S*f S*f SE S*F S*E SE S*E S*E 
S*x*f', 'headerlines',20); 

elseif (type=='P'|type=='p') 

[L photo]=textread(file,'%S*s %S*f S*f S*£ SE S*F S*E S*E S*F S*FE 
Sf', 'headerlines',20); 

elseif (type=='x'|type=='X') 

L photo]=textread(file, 'sf%f','delimiter','\t'); 




















ype mus be either TONEN 3 "ae "PM or Wott 
L=NaN; 
photo=NaN; 
return 
end 





[L2,R2]=textread('am0O1.txt', 'Sf%Sf', 'headerlines',1); 


[a b]=max (photo); 
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R2norm=R2/max (R2); 

photonorm=photo/a; 

plot (L,photonorm, L2,R2norm,L(b),photonorm(b),'rh") 
title('Frequency Response vs. AMO0') 

xlabel ('Wavelength (um) ') 

ylabel('Normalized Photogeneration') 
axkis:().227.3%.0;021)) 

legend('Frequency Response', 'AMO') 


L (b) 























F. TIME.M 


function [T]=time(s) 

This funcion takes an input time, 's', 
Sin seconds and returns time in the format 
Sof hr:min:sec 








%(c) 2004 B.P. Davenport 





hr=floor(s/3600); 
min=floor(s/60-hr* 60); 
sec=floor(s-min*60-hr*3600) ; 
if (hr<10) 
hr=sprintf('O%u',hr); 
else 
hr=sprintf('%u',hr); 
end 
if (min<10) 








min=sprintf('O%u',min); 











min=sprintf('%u',min); 


if (sec<10) 
sec=sprintf('0%u',sec); 
else 
sec=sprintf('%Su',sec); 
end 


T=sprintf('%s:%s:%s',hr,min, sec); 
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